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ABSTRACT 

We present observations and analysis of the host galaxies of 23 heavily dust-obscured gamma-ray 
bursts (GRBs) observed by the Swift satellite between the years 2005-2009, a sample representing all 
GRBs during this period with an unambiguous host-frame extinction of at least Ay > 1 mag. Deep 
observations with Keck, Gemini, VLT, HST, and Spitzer successfully detect the host galaxies and 
establish redshifts for all 23 events, enabling us to provide measurements of the host stellar masses, 
star- formation rates (SFRs), and mean extinctions. Compared to the hosts of unobscured GRBs at 
similar redshifts, we find that the hosts of dust-obscured GRBs are (on average) more massive by about 
an order of magnitude and are also significantly more rapidly star-forming and more dust-obscured. 
However, while the inclusion of this population of dust-obscured hosts shows that GRBs populate all 
types of star- forming galaxies including the most massive, luminous systems at z ~ 2, at redshifts 
below z < 1.5 the overall GRB population continues to show a highly significant aversion away from 
massive galaxies and a corresponding preference for low-mass systems relative to would be expected 
given a purely SFR-selected galaxy sample. This supports the notion that the GRB rate is strongly 
dependent on metallicity, and may suggest that the most massive galaxies in the universe underwent 
an significant transition in their chemical properties approximately 9 Gyr ago. We also conclude that, 
based on the absence of unobscured GRBs in massive galaxies and the absence of obscured GRBs 
in low-mass galaxies, the dust distributions of both the lowest-mass and the highest-mass galaxies 
are relatively homogeneous, while intermediate-mass galaxies (~ 10 9 Mq) have much more diverse 
internal dust properties. 

Subject headings: gamma-ray bursts: general — galaxies: star formation — dust: extinction — ISM: 
structure 



1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs 12 ) represent 
a rare, violent endpoint of stellar evolution. A GRB is 
generated when a newly formed compact object (a neu- 
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12 "Long-duration" GRBs are generally denned as events with 
durations of T90 > 2 s. These are distinguished from short-duration 
GRBs, which typically have T90 < 2s and appear to have a com- 
pletely different origin, most l ikely associated wi th the merger of 
compac t objects (e.g. , INaka r 2007; Bcrgcrl IMTlT ) . Among Swift 
bursts (Gchrcls ct al. 2004) with detected afterglows, observed 
long-duration events outnumber short-duration events by more 
than 10 to 1, and for simplicity, in the remainder of the paper 
we will use the term "GRB" to refer only to the long-duration 
class. 



tron star or black hole 13 ) within a massive star is able 
to briefly power a relativistic jet that pushes through the 
stellar enve l ope and into th e circumstellar medium (e.g., 
lUsovl 119921: IWooslevi [1993ft . Both the radiation associ- 
ated with the ejecta during the explosion (the prompt 
gamma-ray and X-ray emission) and the longer-lived 
multiwavelength afterglow that follows (the emission pro- 
duced by the relativistic shock wave t hat results from 
this explosion; iRees fc Meszaroslll992t iSari et al -111 9 9 8ft 
are extremely luminous, so they are detectable out to 
high cosmolo gical redshifts (including a few at redshift 
z > 8; e.g. iTanvir et al.l 120091: ISalvaterra et al.l 120091: 
ICucchiara et all 12011ft ' even in short observations with 
mid-sized telescopes. 

The association between GRBs and the destruction 
of massive, short-lived stars (e.g., iGalama et all 119981 : 
iHiorth et afll2003t iWooslev fc Hegerll2006ft predicts that 
GRBs should form exclusively in star-forming environ- 
ments. This prediction generally seems to be upheld — 
the hosts of GRBs are ubiquitously young and essen- 
tially always show evidence of recent star formation, both 
in an integrated sense (e.g., ISavaglio et all [20091) and 
at the precise locati on of the GRB within the galaxy 
(jFruchter et alll2006ft . 

A more complicated, still unsettled question is whether 
GRBs form in all star-forming environments — or at 

13 More exotic compact-obj ect constructs such as quark stars 
have also been considered; e.g., Paczynski & Hacnscl (2005). 
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least, if they do so in proportion to the star-formation 
rate (SFR) as one might naively expect given their 
direct association with young stars. While individ- 
ual GRBs are certainly useful probes for studying in- 
dividual star-forming galaxies, a direct, linear associ- 
ation with cosmic star formation would add tremen- 
dous statistical power to studies of the GRB host pop- 
ulation, directly constraining the relative cosmic SFR 
in galaxies of different types (i.e., as a function of 
host mass, luminosity, extinction, morphology, etc.) 
as well as its evolution with redshift, even in galax- 
ies well below the detect i on limit of most flux-limited 
samples (iNataraian et~aLlll997t Hogg fc Fruchterl 119991; 
Diorgovski et al.l 12001b iFvnbo et alj l2001t iBerger et al.l 
2003: IJakobsson et alj|2005l: ITanvir et alll2012D 



Theoretically, there are many reasons to expect that 
this ideal may not be met in reality. Metallicity in partic- 
ular is thought to play an important role in massive-star 
evolution: metals provide opacity to the stellar envelope, 
helping an evolved star to expel its diffuse hydrogen and 
helium layers, which must be removed by some means 
to enable the jet to escape the star and to be consis- 
tent with the obser vation of hydrog en-free Type Ic su- 
pernovae (SNe; see iFilippenkotl 19971 for a discussion of 
SN classification) associated with GRBs 14 . Metals also 
help strip the star of angular mo mentum, a process which 
may inhibit the central engi ne (Wooslcv & MacFadycn 
ll999tlWooslev fc Hegerll2006t ). Given the large variation 
in average metallicity between different galaxies, these 
effects could easily be responsible for large deviations 
between the SFR and GRB rate in different galaxies 
(|Hirschi et al.l I20051) — and metallicity does indeed seem 
to affect the relativ e numbers of different types of SNe 
(jArcavi et al.l l2010h . including the Type Ic broad-lined 
SNe which accompany GRBs. 

Metallicity need not be the only factor: recent evidence 
for variation in the ini ti al ste ll ar mass function (IMF; 
van Dokkum fc Conrovl 120101 : IConrov fc van Dokkunj 



2012f) would also be predicted to produce some vari- 
ation of the GRB rate in relation to other tracers, 
since GRBs are probably produce d only by extremely 
massive stars (at least 20-50 M ; iMazzali et al.1 120031 : 
lOstlin et al.l l2008h . If the distribution of other stellar 
initial properties beyond mass alone (such as rotation or 
binary separation) exhibits similar dependencies on en- 
vironment, these could (in principle) also affect the GRB 

rate relative to that of overall st ar formation. 

While some st u dies (e.g., IJakobsson et al.1 120051: 
Fvnbo et al.l 120081 IChen et ; all [20091 iMannucci et al.1 



20111 : IMichalowski et al. I2012cf) do support the notion 



that the GRB rate is consistent with a model in which 
it is strictly proportional to the overall SFR in some 
situations, the notion that the GRB rate shows sig- 
nificant deviations from the prediction of a uniform 
GRB-to-SFR ratio as a function of environment has re- 
ceived observational support from a number of other 
studies. For instance, the number of GRBs reported 

14 Alternatively, some models enable the star to evolve homo- 
geneously by continuously mixing the outer layers into the core, 
thereby converting the e ntire envelope to heavier elements over 
the course of its lifetime (Yoon & Langcr 2005; Wooslcv & Hcgcr 
2006). Very massive stars may also be able to shed their outer en- 
velopes in m assive eruptions having little to do with the envelope's 
line opacity dSmith fc Owockill2006l) . 



within spiral galaxies appears to be much lower than 
predicted given the amount of total star formation (as 
traced by the rates of Type II SNe) happening in 
these galaxies (|Fruchter et al.l 120061 : IWainwright et all 
120071 ). while the number of GRBs in extremely low- 
mass, low-metalli city systems appears to significantly ex- 
ceed predictions (IStanek et alll2006l; iModiaz et aTll2008l; 
iLevesaue et al.l l2010al) . At higher redshifts, a signifi- 
cant dearth of GRBs within luminous infrared galax- 
ies (LIRGs; Lib, > 10 11 Lq) has also b e en re ported: 
for example, the work of iLe Floc'h et all (|2006l ) found 
only three LIRG hosts from a sample of 16 GRBs ob- 
served, when in reality such galaxies are thought to be 
responsible for about half of all star formation at z > 1. 
Similarly, only a handful of GRBs within ultraluminous 
infrared galaxies (ULIRGs; Lir > 1 12 L(^) or submil- 
limeter galaxies (S MGs) are known (jBerger et al.l 120031 : 
ITanvir et afl2 004). even though these systems contribute 
substantially to, and may domina te, the cosmic SFR 
at higher redshifts (z > 1.5: e.g.. iSmolcic et al.l 120091 : 
IPerez-Gonzalez et alll2005l: IMichalowski et al.ll2010D . " 

However, most of these studies suffered from a signif- 
icant limitation. The ability to search for host galaxies 
is limited by the need to localize the GRB to subarcsec- 
ond precision, and while this is possible at many wave- 
lengths (radio, millimeter, near-infrared [NIR], optical, 
and X-rays are all frequently employed), the majority of 
afterglow positions before the launch of the Swift satel- 
lite were provided optically 15 . For events only observed 
at optical wavelengths, the presence of significant inter- 
stellar extinction within a GRB host galaxy could easily 
conceal the optical afterglow and therefore prevent iden- 
tification of the host. If the dust properties of the GRB 
sightline correlate with those of the host galaxy itself, 
these dust-obscured GRBs could potentially hide an en- 
tire class of host galaxy with properties quite similar to 
those that were largely "missing" from these pre- Swift 
works (massive, luminous, and dusty). 

Dust-obscured GRBs do exist and even appear to 
be fairly common, manifesting themselves as so-called 
"dark" GRBs, events with abnormally faint (and as 
a result usually undetected) optical afterglows 16 . The 
existence of such events has been known almost sinc e 
the beginning of the afterglow era (|Groot et al.l [19981) . 
with dust extinction a l eading candidate for th e ir ex- 
planation (iTavlor et al.l 119981: iDiorgoyski et al.l 120011: 
Klose et al.ll2003HGorosabel et alll2003l : IJakobsson et al l 



2004( 1. More recent work based on deep multi-color 



optical follow-u p observations of l a rge, unbiased sam- 
ples of GRBs (|Cenko et all 12001 IPerlev et al.l l200l 
IGreiner et al.l 120111) has confirmed these earlier suspi- 
cions: approximately 25% of all Swift GRBs are too faint 
to detect even if followed up immediately with a 2 m-class 
ground-based telescope, and while a number of factors 
affect optical detectability, the predominant cause for 
the optical nondetection of bright, well-observed bursts 

15 See, for example, the statistics in th e table complied at 
http://www.mpe.mpg.de/~jcg/grbg emhtmll . 

16 "Darkness" can be defined more quantitatively in various 
ways, and several (conflicting) definitions are employed in the lit- 
erature (see [[2] for an expanded discussion). Here we use the term 
more loosely to denote events with atypically faint optical after- 
glows relative to other wavelengths or to other GRBs at the same 
epoch. 
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is dust extinction, which affects more than 50% of events 
at least modestly and renders perhaps 15-20% of events 
undetectable under normal conditions iPer lev et al.ll200ll 
iGreiner et al.l 1 2 1 ll (A smaller fraction — about ~ 5- 
10% — of all Swift events originate at z > 6; the optical 
undetectability of these bursts is a result of absorption by 
the neutral intergalactic medium.) Studies of the X-ray 
attenuation of optically dark GRBs have reached simi- 
lar conclusions (iFvnbo et al.l2009btlMelandri et al.ll2012t 
IWatson fc Jakobssonll2012l) . It is of obvious interest to 
closely examine the host galaxies of these dust-obscured 
events, both to determine whether they have the poten- 
tial to alter our conclusions about how the GRB rate is 
connected to the cosmic SFR, and to extend the detailed 
analysis permitted by GRB host studies to a wider range 
of environments than those probed only by unobscured 
sightlines. 

Fortunately, the ability to localize events without a 
bright optical afterglow has dramatically improved over 
the past decade. Ground-based follow-up capabilities 
have substantially developed since the pre- Swift era, 
and the early nondetection of an optical afterglow reg- 
ularly motivates deeper follow-up studies from 8 m-class 
telescopes (frequently at NIR wavelengths) or observa- 
tions in unobscured parts of the electromagnetic spec- 
trum (X-ray, submillimeter, or radio) that do success- 
fully detect a counterpart. But even if no additional 
detecti ons are secured. Swift 's onboard X-ray Telescope 
(XRT; iBurrows et al.1 [2005) effectively gu arantees a ~ 
2 position for every long-duration burst (|Butlerl 120071 : 
iGoad et al.ll2007t lEvans et al.ll2009[ ). sufficient to localize 
a host galaxy with reasonable confidence a large major- 
ity of the time. Earlier X-ray cameras on GRB-detecting 
satellites had much poorer angular resolution, and were 
not sufficient for uniquely identifying the host galaxy. 

Studies to date have provided mixed results regard- 
ing the extent to which the host population unveiled 
by dark GRBs actually differs from the host popula- 
tion identified by optically br i ght GRBs. The pre- 
Swift sample of lLe Floc'h et ail ( 20061 ) did include a few 
dark GRBs, only one of which had a sufficiently large 
stellar mass or mid-IR dust emission to be detected 
in Spitzer Space Telescope (Spitzer) observations. The 
darkest GRB host s within the uniform Swift sample of 
ICenko et al.l (|2009t ) did not have markedly unusual opti- 
cal properties disti nguishing them from other GRB hosts 
(| Per lev et al.ll2009T ) , and th e late-time ob s ervatio ns of the 
host of GRB 060923A by iTanvir et all (|2008bf ) identi- 
fied only a faint host with fairly ordinary color. How- 
ever, the hosts of several other individual dark GRBs 
have now shown properties that are actually quite unlike 
those of the supposedly "typical" lo w-mass, \ow-Ay host 
galaxy . Specifically, GR Bs 051022 (ICastro-Tirado et all 
2007bt iRcTeTail [2007ft . 80207 (ISvensson et al.l l2"oia 



Hunt et al.ll201~ll). 08 0325 jHashimoto et al.ll20ltih . and 



080607 dChen et all l20ll have now all been associ- 
ated with quite luminous and massive host galaxies, 
and the host of GRB 020819 h as recently been show n 
to be a high-metallicity spiral (jLevesaue et all l2010bf ). 
While the hosts of optically bright bursts are occasion- 
ally quite luminous in the ultraviolet (UV) and can 
also be moderately massive and occasionally met al rich 
(e.g.. Ilevesaue et all l2010at iKriihler et al.ll2012al) . sub- 
stantial Swift and pie- Swift host surveys and compila- 



tions (|Chen et al.ll2009HSavaglio et alj|2009j ) have turned 
up only a few host galaxies with stellar masses or average 
dust extinctions even approaching those of the dark GRB 
hosts above. Still, given that the study of these indi- 
vidual objects and publication of the resulting discovery 
was in some cases surely motivated by the properties of 
the host itself, it is difficult to determine whether these 
events represent the "typical" dark-burst host or are rare 
exceptions to th e blue-and-fain t rule. 

The study of IKriihler et al.l (|2011h was the first to 
extend these investigations of individual objects to the 
broader population, using a sample of eight dark GRBs 
from the literature and from the Gamma-Ray Burst Op- 
tical Near-IR Detector (GROND), supplemented by New 
Technology Telescope (NTT) and Very large Telescope 
(V1T) observations. Mirroring the results for individ- 
ual dark GRBs, they measure a wide range of properties 
among the hosts in this sample, ranging from relatively 
small, minimally obscured galaxies typical of non-dark 
hosts up to very luminous, massive, and dusty galaxies. 
While their sample size is too small to make strong statis- 
tical statements, these results suggest that massive hosts 
are indeed reasonably common (if not ubiquitous) among 
dark GRBs. Similar results (i.e., a substanti al fraction of 
very r ed galaxies) were seen in the study of iRossi et all 
1)20120 , which targeted 17 GRBs with no detected op- 
tical afterglow. A much larger sample of 69 uniformly 
selected host s obser ved with the V1T was presented by 
iHjorth et aD (|2012f ). While only two-color (R and K) 
photometry of each host is available and deep constraints 
on the presence or absence of an optical afterglow are not 
always available, bursts with no detection of an optical 
afterglow do seem to have substantially redder hosts on 
average (Malesani et al. 2013, in prep.), suggesting a 
trend toward dustier and more massive hosts. 

It is necessary to continue moving from a regime dom- 
inated by studies of individual objects toward statistical 
examination of large samples based on understandable 
selection effects. To this end, over the past several years 
we have been conducting a comprehensive, multiwave- 
lcngth campaign devoted to the observation and charac- 
terization of optically dark GRBs and their hosts. Here 
we present the first results of our campaign — confirming 
that massive, dusty, luminous hosts are in fact typical 
(but not ubiquitous) among the population of "dark" 
bursts, and examining in detail the implications of this 
discovery for the overall population of GRB hosts and 
for the connection between the GRB rate and SFR. 

An outline of this paper is as follows. In $2] we briefly 
describe our selection of the sample. In Jj3]we present our 
observations of the hosts with Keck, Gemini, Spitzer, and 
the Hubble Space Telescope (HST). Section U describes 
our method for estimation of stellar masses, SFRs, and 
minimum bolometric luminosities from these measure- 
ments and other observations in the literature. In Sj5] 
we outline each burst individually; we describe the char- 
acteristics motivating its inclusion in the sample, iden- 
tify its host galaxy and discuss observations thereof, and 
summarize the host's properties. We examine in Sj5]thc 
properties of the entire sample as a population in com- 
parison to other, previously published samples of GRBs 
and other populations of high-z galaxies. In fj7] we dis- 
cuss the implication of our results for the origins of dark 
bursts, the distribution of dust in high-z galaxies, and 
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Fig. 1. — Histogram of measured or limiting afterglow extinc- 
tion columns (in terms of rest-frame Ay) of several samples of 
Swift GRBs, including events presented in this paper. The un- 
filled yellow-and-black line shows an estimate of the intrinsic, un- 
biased distribu tion based on a ne arly complete sample of events 
combined from lCenko et al.1 lf200l ) and lGreiner et all if^oTH ). The 
hatched blue histogram shows measured values of Ay for a s ample 
of well-observed , optic ally bright events from [Kami et al. (2006) 
and Kami et al. (2010), rescaled to match the < 0.25 mag bin for 
the complete sample; the ~ 20% dustiest events are systematically 
missed in these optically selected studies including effectively all 
GRBs with Ay > 2 mag. The solid red histogram shows the Ay 
measurements (or lower limits) for the 23 events analyzed in this 
sample. 

the ability of GRBs to serve as unbiased probes of star 
formation. Our conclusions are summarized in fj8] 

2. SAMPLE SELECTION 
2.1. Motivation 

Our primary goals in this paper are to characterize the 
GRB host population missed in earlier work due to dust 
extinction of the optical afterglow and to examine the 
implications for inclusion of this missing population on 
larger questions of the coupling between the GRB rate 
and SFR, and on general questions relating to GRB hosts 
and high-redshift galaxies. Approaching this task is not 
straightforward, since optical coverage of GRBs is spotty 
due to the unpredictable times and locations of these 
events — GRBs often occur too close to the Sun or Moon 
to observe, and ground-based work may by stymied by 
weather and other conditions. Many optically bright 
bursts are missed due to these prosaic reasons, and with- 
out careful attention to selection (for example, had we 
merely chosen a sample of events with no reported after- 
glow) our sample may be "polluted" by ordinary events 
and lead us to incorrect conclusions about the types of 
events that are systematically missed. 

In addition, many factors other than dust extinction 
can influence the optical brightness (and therefore de- 
tectability) of a particular GRB at a given time: the 
burst's energetics, its circumstellar density, the temporal 
evolution of its light curve, as well as a variety of micro- 
physical parameters and of course the burst's distance 



and redshift (e.g.. ISari et al.l 1998 


: iGroot et al.l 119981: 


Diorgovski et al. 2001; Fvnbo et al. 


2001: lLazzati et al.l 


20021 iNvsewander et al.1 12009D. These causes are all 



physically distinct, and most of these factors proba- 
bly have little to do with the large-scale properties 
of the host galaxy the burst inhabits. On the other 
hand, given the strong correlations observed between the 
mean obscuration of a galaxy's stellar population and 
its fundamental properties such as mas s and SFR (e.g., 
iMeurer et ail 119991 : IShaplev et aTll2001l ). there is signif- 
icant reason to expect that dust-obscured GRBs may 
indeed reveal a different host population. Dust obscu- 
ration also is the predominant cause of optical nondc- 
tcction within s amples where GRB s are well observed 
at early times (ICenko et al.l 120091: iPerlev et all 120091: 
iGreiner et al.ll2011l ). especially among luminous and well- 
observed populations. 

For these reasons, in this study we focus only on GRBs 
whose afterglows have been heavily absorbed by dust in 
the interstellar medium (ISM) of their host galaxy. We 
will generally refer to this class as dust-obscured GRBs to 
distinguish them from the more general class of optically 
"dark" bursts, which can result from numerous causes 
(or combinations of causes), although we will occasion- 
ally continue refer to dust-obscured events using simply 
"dark" as a shorthand. 

2.2. Selection Criteria and Implementation 

Since we desire a quasi-complete sample of all known 
dust-obscured events during a given time period to avoid 
biasing our host sample towards galaxies that are partic- 
ularly bright or have other conspicuous properties that 
may have attracted our (or others') attention first, we re- 
strict our search at the outset to events observed by Swift 
that occurred during the five-year period of 2005-2009. 

Within this temporal window, our primary condition 
for the inclusion of a given GRB in our sample is direct 
evidence (from observations of the X-ray, optical, and 
NIR afterglow) for extinction of at least Ay = 1 mag in 
the host frame 17 . We select this threshold on the basis of 
th e fact that that very few optically-bright GRBs exc eed 
it (|Kann et al.ll2fM WM . iSchadv et al1l2007l HoT3) 18 , 
yet events of higher Ay do represent a significant contri- 
bution to the overall GRB rate (~ 15%; Figure 1). GRBs 
with Ay > 1 mag therefore represent a population that 
is intrinsically common yet highly underrepresented in 
previous host work. This magnitude threshold is also 
achievable in practice: as long as a burst is followed up 
rapidly with a moderate-size telescope (or within the first 
day by a large-aperture telescope) it is usually possible 
to determine whether it has Ay < 1 mag or Ay > 1 mag. 

Since for heavily obscured events the optical or NIR 
afterglow is often not detected at all, making this de- 
termination requires a few assumptions on th e intrin- 
sic spectrum. Following Jakobs son et al.l (|2004| ). we as- 
sume that base d on simple synchrotron models (e.g., 
I Sari et all 119981 ) a fading GRB afterglow must have 
an intrinsic spectral index (defined as F v oc v~P) of 

17 Throughout this section we remove G alactic foreground ex - 
tinction at the outset using the dust maps of Schlcecl et al. ( 1998) 
and neglect intergalactic extinction. 

18 Although 1 mag of attenuation is a relatively small amount 
and GRBs are very luminous, since the typical GRBs lies at z = 1-3 
this corresponds to much greater attenuation in the observer frame. 
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Pox > 0.5 in all circumstances, or (more stringently) 
an intrinsic /?ox > Px.min — 0-5 if late-time XRT ob- 
servations ()Butler fc Kocevskil 120071 ) indicate an X-ray 
spectral s lope steeper than 0x = 1 at 95% confidence 
(following Ivan der Horst et~aIll2009L but we use the 95% 
confidence lower-limit on Bx rather than the best-fit 
value). 

For events where there is no detection of an opti- 
cal/NIR counterpart (or there is one detection in only a 
single NIR filter, typically the K band), we simply take 
this baseline minimum intrinsic flux and determine the 
minimum extinction needed beyond this to explain the 
optical/NIR nondetections. The Swift XR T X-ray flux 
is determin ed from the Swift data pages (|Evans et al.l 
120071 12009T ) , converted to a 1 keV flux density calculated 
using the average spectral index and absorption correc- 
tion, and interpolated to the exact time of observation 
based on a power-law fit to neighboring data points. An 
extinction curve like that of the Small Magellanic Cloud 
(SMC) 19 (fixed either to the known redshift or a fiducial 
z = 2 if the redshift is unknown) is then applied to deter- 
mine how much dust is necessary to interpret the upper 
limit or measurement. 

Even if analysis of the optical/NIR points individually 
versus the X-ray flux does not unambiguously establish a 
large Ay, a separate constraint can be established based 
on the optical or NIR color — an extremely red afterglow 
can potentially provide much stronger (and physically 
definitive) constraints on extinction than even a quite 
deep nondetection. If the afterglow is detected in at least 
one filter, all UV/optical/IR points and limits are scaled 
to the same epoch (using the optical light curve if there 
are enough observations in the same filter to establish it; 
otherwise by the optical decay rate in the range < a < 
2 that produces the bluest spectral energy distribution 
(SED)), and we again try to determine the minimum Ay 
that can simultaneously explain these observations (and 
the X-ray observations) within the synchrotron model. 
Both an unbroken intrinsic optical-to- X-ray SED and an 
SED with a break between the bands (0o = 3x — 0.5) 
are attempted. 

In principle, we must take into account the possibility 
of Lyman-cc or Lyman-break absorption from the inter- 
galactic medium in this analysis, since these can also 
significantly attenuate the flux in blue or UV filters or at 
high redshift. At most redshifts typical of GRBs this is 
not a concern, since in practice the redder filters (R band 
and redward) always provide the strongest constraints on 
extinction anyway, and the R band is not significantly 
suppressed at z < 5. However, we do have to consider 
the possibly that a given GRB is at z > 5 as an alter- 
native hypothesis; higher- z GRBs represent a small but 
nontrivial fraction of the GRB sample (about 5-10%; 

19 The SMC curve is employed because it produces the largest 
rest-frame UV (observer-frame optical for typical Swift GRB red- 
shifts) extinction for a given Ay among the Local Group extinction 
curves, and therefore is the most appropriate choice for evaluating 
the minimum Ay. Of course, even steeper extinction curves, such 
as those im plied from populations of Type la SNe and perhaps 
some GRBs IPoznanski et al.ll2009l ; IButler et a"I]|2006l ; IZafar et all 
120111) , could produce (slightly) lower values of Ay for a given Atjv , 
so in this sense our Ay limits are not strict. Nevertheless, for con- 
sistency and simplicity we use the SMC curve in this work in cases 
where the extinction curve cannot be inferred directly from multi- 
color detections of the afterglow. 



e.E.JPerlev et al.l2009 
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120111: Uakobsson et al. 


20121). We can exclude the high- 



redshift hypothesis in any of the following ways. 



1. From detection of the transient in any optical fil- 
ter bluer than the filter (s) used to determine the 
presence of significant Ay. (The combined effects 
of Lyman-break and Lyman-a absorption produce 
a sharp cutoff at z > 5 and partial attenuation is 
only possible in one broad-band filter at a time.) 

2. From a measured /3ox < 0.5 in the K or H bands, 
which are not absorbed by neutral gas except at 
extremely high redshifts (z > 11). 

3. From detection of significant X-ray absorption in 
excess of the expe cted Galactic value. Following 
iGrupe et al.l ()2007f ). we determine the Ah excess 
(fit from the XRT data at z = with the known 
Galactic column subsequently subtracted) to rule 
out very high redshifts. Spec ifically, we employ 
the au tomatic online tables of IButler fc Kocevskil 
([200l 20 and apply a minimum threshold of Ah > 
10 21 cm~ 2 to the lower limit on Ah in that work. 
(The value of 10 21 cm -2 correspon ds to a limit 
of Zma x = 5 in the formula given by IGrupe et al.l 
(|2007l ). While this criterion was developed based 
only on events through 2007, we empirically ver- 
ified that it still comfortably excludes all known 
z > 5 events among Swift GRBs through to the 
present time.) 

In all cases, the lower- redshift association was indepen- 
dently verified by the detection of a host galaxy under- 
lying the afterglow position and direct measurement of 
its redshift spectroscopically or photometrically. (How- 
ever, host detection and redshift measurement was not a 
criterion for inclusion.) 

The implementation of this procedure in practice will 
be described in detail in a separate work (Perley et al. 
2013, in prep.) In brief, we downloaded all photomet- 
ric observations (including upper limits) from the GCN 
Circulars in 2005-2009 as we ll as the onl i ne library o f 
XRT light curves provided bv lEvans etldl (|2007l l200ll . 
and automatically calculated Box and Ay^ m { n for each 
point, as well as Ay^ t if possible. Events close to satis- 
fying the Ay > 1 mag criterion were followed up in more 
detail with observations from published or unpublished 
sources as necessary (where possible) in order to verify 
this association and provide the deepest possible limit. 

A few targets were excluded (specifically GRBs 050716, 
060923C, and 080229A) due to the presence of a bright 
Galactic star within 2" of the GRB afterglow posi- 
tion. We also excluded GRB 070412, which is in the 
outer halo of a bright foreground galaxy. Moreover, we 
elected not to include GRBs 060807, 080605, 080805, and 
090926B, for which certain estimates indicate a best- fit 
value slightly above Ay « 1.0 mag but with a range of 
uncer tainty (or alternate solutions) permit ting a lower 
value (jKruhler et al.ll201lHZafar et aLll2012ft . After these 
subtractions, the final Ay > 1 mag sample contains a to- 
tal of 23 GRBs as summarized in Table [T] In brief, the 
sample has the following properties. 

20 http: //butler . lab . asu. edu/Swif t/xrt_spec_table .html 
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• All but one event (GRB 051022) were detected 
by Swift. While we do not restrict our search to 
Swift bursts, our procedure requires XRT follow- 
up observations nearly simultaneous with the opti- 
cal measurements in order to evaluate the spectral 
index. 

• Eight events were detected in multiple NIR/optical 
filters and can provide a direct estimate (that is, 
well-constrained lower and upper limits) on the 
rest-frame extinction. Six were detected in only one 
filter (always only the K band), and while an up- 
per limit could be placed in principle by assuming a 
minimum redshift and maximum Pox , it would be 
very large (typically ~ 10-30 mag) and not partic- 
ularly useful. The remaining ten were not detected 
in any NIR/optical filter and have only lower limits 
on Ay. 

• All but four events (GRBs 060923A, 070802, 
071021, and 080325) have detection of significant 
./Vh excess. Among the exceptions, GRB 070802 
has a spectroscopic afterglow redshift and the re- 
maining three had if-band afterglow detections be- 
low /?ox = 0.5, so they are also extinguished (as 
opposed to being at high redshift). Indeed, addi- 
tional observations confirmed the presence of a host 
galaxy underlying all four positions. 

The prompt-emission and afterglow properties of our 
sample are also summarized in Table [T] and in Figure 

m 

2.3. Impacts of Possible Biases in Sample Selection 

Note that our sample is necessarily not complete, con- 
stituting 23 events out of ^400 long-duration GRBs lo- 
calized by the XRT during the five-year window (among 
which ~ 40-80 probably had a "true" Ay > 1 mag). The 
majority of dust-obscured bursts are inevitably missed, 
because the rapid or deep follow-up observations neces- 
sary to identify them conclusively are usually not con- 
ducted. In addition, the inclusion or exclusion of a few 
events close to the Ay ~ 1 mag threshold could be de- 
bated (i.e., varying assumptions about /?ox, the choice 
of extinction law, the X-ray fitting procedure, etc. would 
alter a few events in the sample) . Neither of these points 
is problematic for our sample, which aims only to gather 
a set of events that is representative of the population 
missed by optical-afterglow searches without introducing 
a dependence on the characteristics of our host galaxies. 
(Such a bias would, for example, be present if we had re- 
quired a reported host redshift, which would necessarily 
disfavor faint galaxies.) 

As the vast majority of afterglow observations were 
carried out (and reported in the GCN Circulars) before 
the host galaxy itself was identified, our procedure above 
should be almost completely independent of the proper- 
ties of the host. A weak bias in favor of brighter hosts 
could in principle be present only as a result of the fact 
that in a some cases host-redshift measurement preceded 
selection, and in principle some events would not have 
made the Ay > 1 mag cut if the fiducial z w 2 was as- 
sumed at the outset instead of the actual redshift. In 
practice, only GRBs 060202 and 090417B would have 



failed our cut at z « 2 while passing it at the host- 
measured redshift. 

The sample is clearly not unbiased in terms of the in- 
trinsic afterglow properties — in particular, events with 
more luminous X-ray afterglows will naturally be favored 
as optical or NIR follow-up observations to a given depth 
are increasingly likely to be constraining if the X-ray 
afterglow flux is higher. For a given X-ray luminos- 
ity, events with relatively flat intrinsic spectral indices 
(j8 ~ 0.5) will be weakly favored if the extinction is very 
large since relatively less extinction is needed to sup- 
press the afterglow, but events with steep intrinsic in- 
dices (j3 « 1) can be favored if the extinction is more 
modest, since optical detection of a reddened afterglow 
becomes feasible. However, as these properties are intrin- 
sic to the GRB itself (the spectral index is also sensitive 
to the immediate circumburst density), we do not ex- 
pect strong correlations with the large-scale host-galaxy 
environment. We therefore anticipate that our sample 
should be reasonably representative of the hosts of "all" 
Ay > 1 mag bursts, including those with fainter after- 
glows. 

We also expect some biases with redshift, since it is 
easier to place constraining limits on the afterglow of an 
event that is nearby relative to an event with similar 
luminosity that is far away 21 . Because rest-frame opti- 
cal/NIR observations correspond to bluer light at higher 
redshifts where even relatively small dust columns will 
absorb a large amount of light, very large extinction val- 
ues also become harder to recognize for this reason (es- 
pecially at z > 4). For these reasons we certainly cannot 
expect our sample to provide a representative redshift 
distribution of all dust-obscured GRBs. Nevertheless, as 
long as comparisons are restricted to objects at similar 
redshifts, this possible bias should not affect any conclu- 
sions drawn by the set of hosts probed by our sample. 

3. OBSERVATIONS 

We observed the fields of all 23 targets satisfying the 
above criteria using a variety of resources from both the 
ground and space. In the following sections we briefly 
summarize these observations, as well as the reduction, 
calibration, and analysis of the data. 

3.1. Keck/LRIS 

The Low Re solution Imaging Spectrometer (LRIS; 
lOke et ail 119951) on the Keck I telescope is an optical 
imager and spectrograph equipped with both blue- and 
red-optimized cameras split by a dichroic. Imaging ob- 
servations of the galaxies were acquired mostly as part of 
our multi-year GRB host follow-up campaign and were 
reduced via standard procedures (Perley et al. 2013, in 
prep.). Photometric calibration was perfo rmed relative 
to th e Sloan Digital Sky Survey (SDSS; lAihara et al.l 
120H in cases where SDSS covered th e field in ques - 
tion; otherwise, we used observations of lLand olt (2009) 
standards (on photometric nights) or our own secondary 
standards obtained with either the 1 m Nickel telescope 
at Lick Observatory or the 60-inch telescope at Palomar 

21 On the other hand, lower-luminosity GRBs may not be de- 
tected by the satellite in the first place at greater distances, and an 
afterglow will stay bright for longer due to time dilation, at least 
partially compensating for this bias. 
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Fig. 2. — Equivalent coeval afterglow SEDs for the sample of 23 Swift-observed GRBs whose host galaxies represent the subject of 
this paper. The gray line shows the intrinsic afterglow flux corresponding to the model with minimum Ay that is consistent with the 
observations; the black curve indicates the observable flux for this model including extinction and X-ray absorption. The dashed line 
indicates /3qx = 0.5 unless the minimum-flux model itself assumed this value. SEDs are normalized to the observed (absorption-corrected) 
X-ray flux at 1 keV. 

A small number of host galaxies were also observed 
with LRIS in long-slit spectroscopic mode. These data 
will not be discussed in detail here, except as they pertain 
to establishing or constraining the host-galaxy redshift. 
A more comprehensive spectroscopic study of dark GRB 
hosts will be presented in future work. 

3.2. Keck/NIRC 

We employed the N ear-Infrared Camera (NIRC; 
iMatthews fc Soiferlll99l on Keck I during a one-night 
classical run on 2009 May 31 (UT is used throughout 
this paper) to observe several objects. These images 
were reduced using standard NIR techniques within a 



Observatory (|Cenko et alj |2006b) . Magnitudes 22 are de- 
termined via aperture photometry, using a custom wrap- 
per around the aper task of the IDL Astronomy User's 
Library 23 . In cases where the host is seen to be resolved 
in any filter we use a consistent choice of aperture for 
all filters, with exceptions for images having particularly 
bad seeing (> 1.0") where a larger aperture is employed. 
Photometry (from LRIS and all other instruments, be- 
low) is presented in Table |3] 

22 Except where specified, when reporting apparent magnitudes 
or colors, we use the Vega sy stem (for non-SPSS filters) or the 
SDSS system (for SDSS filters; iFukiigita et al.|[T99l) . 

2,! http: //idlastro . gsf c .nasa. gov/ . 
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Fig. 3. — Broad-band photometry of our sample of dark GRB host galaxies fitted by our population-synthesis SED modeling procedure. 
Red points are photometry (2a uncertainty error bars), with the black curve indicating the best fit and the gray envelope the lcr model 
uncertainty region. Gray triangles indicate upper limits (2cr). 

custom Python pipeline. Seeing conditions were good 
and the night was photometric, so we used observations 
of the standard stars FS23/M3-193, FS33/GD153, and 
FS29/G93-48 to establish the photometric calibration for 
most fields. We checked these for consistency using stars 
present in our science fields with magn itudes from the 
Two- Micron All Sky Survey (2MASS; ISkrutskie et all 
I2006T) . or with P AIRITEL (the Pet ers Automated In- 
frared Telescope; iBloom et al.l I2006D calibration obser- 
vations which we used to extend the 2MASS calibration 
down to fainter stars, when possible. PAIRITEL employs 
the same camera, telescope, and filters as the northern 
2MASS survey. 



Deep ima ging observations fro m the Near-Infrared Im- 
ager (NIRI: iHodapp et al.l 120031 on Gemini-North were 
used to constrain the properties of host galaxies in the 
sample. Data came both from our target-of-opportunity 
(ToO) program (in cases where no transient behavior was 
observed at early times, or when a late-time image was 
obtained to confirm suspected early variability) and from 
two classical nights in 2010 obtained as part of the Keck- 
Gemini exchange program. Data from both runs were 
reduced using the NIRI reduction utilities in the Gemini 
IRAF package, automated using a custom Python script. 
We calibrated the observations using 2MASS standards 
in the field or against our PAIRITEL calibrations. 



3.3. Gemini/NIRI 



3.4. Hubble Space Telescope 
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We have obtained HST images of several of our targets. 
Our typical strategy (employed in GO programs 11343, 
11840, 12378, and 12764; PI Levan) was to obtain ob- 
servations in one optical and one NIR band, for which 
we chose F606W (broad V to R) and F160W (broad H) 
to provide a good combination of sensitivity and wave- 
length range. Sometimes we acquire two orbits of ob- 
servations (one per filter) to obtain these data, in which 
case we use ACS and WFC3 for the optical and NIR, 
respectively, with a standard 4-point box dither patten. 
In other cases we use a single WFC3 orbit, split between 
F606W and F160W, with a 3-point line dither in each 
filter. The resulting data were combined within mul- 
tidrizzle, with pixfrac=l and the scale left at the na- 
tive value for these moderately dithered images. We also 
supplement these images with data taken from a separate 
Snapshot program (using WFC3-IR in just the F160W 
band; GO-12307) for GRB 070521, plus some imaging 
from our ongoing Cycle 20 program (using WFC3-IR in 
a variety of filters; GO-12949) and archival observations 
from several earlier programs. Images of GRBs 080207 
and 080607 were taken from previous studies of these 
objects (see lChen et~aIll201Ct liSvensson et al.ll2012t ). 

We obtain magnitudes by aperture photometry vi a 
mag-auto within SExtractor ( Bertin fc Arnoutsl [1996): 
results are presented in Table El In cases of observations 
in two filters we use the same physical aperture for con- 
sistency, and in most cases we employ a similar aperture 
as that used in ground-based imaging as well. 

3.5. Spitzer Space Telescope 

Observations of many of our targets were carried out 
with the InfraR ed Array Camera (IRAC; iFazio et al.l 
12004 ) on Spitzer ([Werner et al.ll2004D during Warm Mis- 
sion Cycle 7. In all cases we employed dithered 100 s ob- 
servations of the field, usually totaling 1500 s (15 dither 
positions) per filter, but occasionally more or less de- 
pending on the anticipated magnitude from previous 
ground-based observations of the target. Targets were 
observed in Channels 1 and 2 (3.6 and 4.5 /mi, respec- 
tively.) GRBs 060923 A and GRB 080207 were previously 
observed during the cold mission, in each case in all four 
IRAC filters. 

We downloaded the PBCD images from the Spitzer 
archive to disk and identified the host galaxy by ref- 
erence to our ground-based or HST observations. For 
several fields, the host galaxy is blended with one or 
more nearby sources in the field; the diameter of the 
point-spread function (PSF) of IRAC is ~ 2". In these 
cases, we used the GALFIT package ([Peng et all 120021 ) 
to subtract the nearby objects based on a model of the 
PSF measured from a bright, isolated star elsewhere in 
the IRAC image. Photometry of the host was then per- 
formed using IRAF, employing an aperture radius of 2 
native pixels (4 resampled pixels in the PBCD observa- 
tions, or 2.4") with a sky annulus of inner and outer 
radii 8 and 12 pixels, and calibrated via the zeropoint 
values in the Spitzer IRAC handbook. In all cases, the 
angular size of the host galaxy was sufficiently small that 
aperture effects are not significant. 

3.6. Very Large Telescope 

Many hosts in this sample were observed with the 
Very Large Telescope (VLT). Most of these data were 



previously published; in particular, we rely heavily on 
the FORS2 i?-band and ISAAC i f -band photometry 
products from the TOUGH survey ([Hiorth et all 120121) 
and several X - shooter redshifts reported in the work of 
iKriihler et all (j2012bl) . We also use our independent 
rereduction of the X-s hooter observations reported by 
ISalvaterra et al.l (|2012l ). retrieved from the ESO archive. 

3.7. Previous Work 

We acquired additional photometry from a variety 
of other literature sources, usually burst-specific papers 
including photometry of the host galaxy under study. 
These sources are cited where appropriate. 

3.8. Host Identification 

Historically, the most significant challenge affecting the 
identification and characterization of dark GRB hosts 
has been the difficulty in localizing the afterglows with 
sufficient accuracy (~ 2-3") to uniquely identify the host. 
As we have mentioned, Swift now provides positions 
accurate to this level for essentially all (long-duration) 
GRBs observed with the XRT, based o n th e optical/UV 
registr ation techniques of lButlerl (|2007l ) and lEvans et al.l 



reg: 
(20 



20091 ). 

Nevertheless, XRT positions are never better than 
~ 1" accuracy and occasionally can leave some ambiguity 
about the host identification. In most cases, the after- 
glow was also detected at optical or NIR wavelengths, 
providing (in principle) a position accurate to 0.5" or 
better. These positions were obtained either by our own 
observations, from published work, or from the GCN Cir- 
culars. (In cases where the accuracy of a reported po- 
sition in the GCN Circulars was in any doubt, we reob- 
tained the afterglow image and recalculated the position 
internally.) 

We have also been acquiring ToO observations of opti- 
cally undetected GRB afterglows using the Chandra X- 
ray Observatory (CXO). The ToO observations of ~ 15 ks 
per exposure were typically conducted with ACIS-S in 
imaging mode within a few days of the occurrence of the 
GRB. Coordinates from these observations (placed in an 
absolute frame using the 2MASS catalog) are presented 
in Table[2j along with our optical/NIR or XRT positions 
where appropriate. 

The host galaxy was identified in the standard way by 
finding the brightest source in the image that is consis- 
tent with the afterglow position. In nearly all cases this 
is almost completely unambiguous: a single, moderately 
bright, well-detected object is seen directly at the after- 
glow position. The significance of the association (follow- 
ing the standard P c hance metric for estimating the proba- 
bility that a circle of a given size placed randomly in the 
sky encloses a galaxy b righter than a certain flux level; 
e.g., iBloom et al.l 120021 ) depends on the filter chosen to 
make the comparison — we adopt the K band for our tar- 
gets since it does not penalize very red, dust-attenuated 
galaxies, although the results would not be qualitatively 
different if the more commonly employed R band were 
chosen. (In cases where K photometry was not avail- 
able, magnitudes are calculated by interpolation from 
other filters.) In every case, the probability of false as- 
sociation is low (always P c hanco < 0.1, and in most cases 
-Pchancc < 0.01). Because of the large size of the sample, 
the cumulative probability of having a small number of 
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Fig. 4. — Optical mosaic, showing a 10"xl0" cutout of an image chosen from our optical imaging of eac h source (ge nera lly, the filter 
showing the clearest detection was chosen). Blue circles indicate X-ray positions from the Swift XRT (from Butler 2007 and Evan s et al.l 
2009) or from CXO, and red circles indicate a position at a longer wavelength (usually NIR, but in a few cases radio or optical). The host 
itself is centered in each tile and encircled with a dotted green line. NOTE: Figures 4-6 have been compressed to lower resolution due to 
arXiv file size limitations. 



misidentifications present is significant despite these low 
individual probabilities: the probability of at least one 
misidentified host is 39% and the probability of two or 
more misidentified hosts is 8%. However, as our study 
focuses on the aggregate properties of the sample and 
is not particularly concerned with the properties of any 
individual case, we do not expect this to be a significant 
limitation. 

4. MODELING 

A primary goal of this work is to infer the proper- 
ties of the host galaxy (in particular the UV-based SFR, 



the average extinction Ay, and the total stellar mass 
M„) for direct comparison to previous samples. While 
these parameters can be crudely estimated from individ- 
ual photometric points independently using various pre- 
scriptions (e.g., using the optical slope /3 to measure Ay 
and correcting the rest-frame UV luminosity to measure 
the SFR; or using the rest-frame X-band luminosity as 
an indicator of stellar mass), given the complex interplay 
between the parameters and the large variation of red- 
shifts and photometric completeness across our sample, 
we instead adopt a more general population-synthcsis- 
based SED-fitting procedure to fit all parameters to all 
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Fig. 5. — Near-infrared mosaic of our GRB host targets. Image sizes and overplotted information are the same as in Figure [4] Blank 
tiles indicate events without our own NIR imaging (NiR photometry was taken from the literature for these events.) 



data points simultaneously. 

The fitting of the host SEDs is implemented with 
our own software wri tten in IDL using the p opulation- 
synthesis libraries of iBruzual fc Chariot] ()2003l ) (here- 
after BC03), assuming a IChabrierl ()2003D IMF and the 
Padova (1994) stellar evolution tracks. Taking as in- 
puts the galaxy's total stellar mass, average metallicity, 
current SFR, and an analytic expression for the time- 
dependent star-formation history, the code combines the 
ap propriate BC03 templ ates and extinguishes them using 
the lCalzetti et all (|1994l ) attenuation law to produce the 
galaxy-integrated spectral luminosity distribution. Ma- 
jor UV/optical nebular emission lines are included, cal- 
culated from the SFR and metallicity following the pre- 



scrip tions oflKennicuttJ (1 19981) , iKewlev k Dopital (|2002ft , 
and IKewlev et al J (|2004h . since in some cases these 
can contribute substantial flux even to broad-band fil- 
ters. These are then redshifted appropriately to produce 
the observed SED, and an intergalactic neutral-gas ab- 
sorption prescription (a semiempirical sightline-averaged 
Lyman-a forest model as a fu nction of wavelength and 
redshift following Madau 1995) is also applied. 

The Calzetti law — a featureless curve similar to a 
power law — is used to model the host-galaxy SED even 
in cases where direct line-of-sight measurements through 
the host galaxy provided by the GRB afterglow favor 
a different extinction law, for three reasons. First, the 
Calzetti curve is the standard one employed in most cur- 
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Fig. 6. — IRAC mosaic showing 3.6 fim imaging with Spitzer. Despite our relatively short integrations (typically 1500 s), we detect 
every object in our sample, indicating a significantly higher average stellar mass than in previous, optically selected GRB host samples. 
In a few cases the host is blended with a nearby object; the host flux is recovered by fitting and removing the flux of any possible nearby 
contaminants before performing aperture photometry. 



rent high-z galaxy work (where it has generally been 
successful at modeling the UV SEDs of distant galax- 
ies in current sur veys; e.g., ICalzetti &; Heckmanl Il999t 
iReddv e~ al. 2010), and this convention facilitates direct 
comparison of the Ay and SFR to other work. Second, 
the Calzetti law is explicitly an attenuation law rather 
than an extinction law. An extinction law refers to a 
point source shining through a single sightline (with a 
single optical depth and in which photons are lost to 
scattering), whereas an attenuation law refers to the ra- 
tio of total emitted energy to total energy escaping (the 
exact situation desired here, in which stars can exist at 
varying optical depths in the cloud and in which scatter- 



ing photons are not lost). Third, in cases in our sample 
where sufficient rest-frame UV data were available to at- 
tempt multiple extinction laws, the Calzetti extinction 
law was favored (e.g., there was no evidence of a 2175 A 
absorption feature or strong UV curvature in the host 
SED). 

The SED is then converted to an observable spec- 
tral flux distribution (SFD) using a basic correction 
for distance (we employ Oa = 0.7, Qm = 0.3, Ho = 
70kms~ 1 Mpc^ 1 ) and redshift. Finally, synthetic pho- 
tometry is performed (using a library of standard filter 
transmission curves) to produce the predicted fluxes in 
various standard filters as observed at a given redshift. 
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The broad-band observations can then be fitted against 
our model using the Levenberg-Marquardt least-squares 
fitting procedure as implemen ted in the mpf it package 
l|Morelll978t lMarkwardtll2009D . Upper limits are treated 
as a flux measurement of zero with an Na uncertainty 
equal to the uncertainty level of the limit. (There are 
only a few constraining upper limits in the sample.) 

Our set of broad-band observations is generally too 
limited to actually constrain the star-formation history 
or metallicity independent of other parameters, even for 
the best-studied galaxies in our sample. To help stan- 
dardize our results against each other and against other 
studies, we apply a few simplifying assumptions. The 
stellar metallicity for all galaxies is assumed to be 0.5 Z© 
(Z = 0.01), with two exceptions. For the host of GRB 
061222A, the UV spectral index is too blue for any solar- 
metallicity template and we instead use 0.1 Z Q . For the 
host of GRB 051022, a v alue close to solar ha s been mea- 
sured spectroscopically ([Graham et al.lf2009h and so we 
assume 1.0 Z©. (In general, however, the assumed stellar 
metallicity does not make a significant difference on our 
results relative to other sources of uncertainty.) Second, 
we apply a simple prescription for the star-formation his- 
tory as follows. The system is first fit with a constant 
star- formation history. If the effective age (M»/SFR) 
converges to a value in excess of 1 Gyr or if the fit is 
poor, we allow the young and old populations to be fit 
separately by permitting a step-function change of the 
SFR at 5 x 10 7 yr before present. The age of the older 
population (i max ) is a free parameter in the fit in this 
model, but is constrained to be less than the age of the 
Universe at that redshift. (In one case, GRB 080325, the 
age was forced to this maximum limit, and for this case 
only we modified the star-formation model to be expo- 
nentially decaying with time for the old population.) If 
the effective age converges to less than 1 Gyr, we retain 
the assumption of a constant star-formation history (we 
are not sensitive to variations) but require a minimum 
value for the effective age of > 2 x 10 7 yr. 

To estimate the statistical uncertainties in each pa- 
rameter, we run a series of 100 Monte-Carlo trials for 
each galaxy in which each flux measurement is modified 
by a random amount (following a Gaussian distribution 
with a equal to its observational uncertainty) . Reported 
uncertainties in the output parameters correspond to the 
15th-to-85th percentile (that is, la) of the results. These 
results are given in Table |U Note that these values do 
not include the systematic uncertainties due to assump- 
tions about the extinction, st ar-formation history, and 
IMF, which can be significant (|Michalowski et ai1l2012al ) 
but should affect all galaxies (and comparison samples) 
in a similar way. Our primary interest is in the aggregate 
properties of the various samples in comparison to each 
other (rather than in the absolute values of individual 
objects), and we do not expect large systematics that 
would endanger these comparisons. 

5. GRBS AND HOST GALAXIES 

5.1. GRB 050915A 

GRB 050915A was a relatively faint and fast-fading 
event with minimal late-time follow-up observations. 
However, it was observed rapidly by a number of tele- 
scopes including (critically) the roboticizcd Palomar 60- 
inch telescope (P60) and PAIRITEL. P60 did not de- 



tect the transient in any exposure (|Cenko et all [20091. 
PAIRITEL reported a weak detection of a transien t 
source in an early f/-band stack (Blo om fc Alatalofe OOS) , 
although no detection was reported in simultaneous J or 
K s imaging. We reanalyzed these frames and performed 
photometry on all three filters by forcing an aperture at 
the location of the .ff-band detection and do detect (at 
about 4-5<r) a weak source in the other filters also. The 
NIR colors and optical nondetection are consistent with 
a dust-reddened afterglow with an inferred extinction of 
Ay ~ 1.4 mag. 

Previous i magin g of this field was reported by 
iPerlev et al.l ([2009) , including the identification of a host 
galaxy offset by 0.85" from the transient location. The 
photometry presented in Table [3] and used in mod- 
eling has been reanalyzed since that time using our 
Lick Nickel 1 m field calibration. Combining these ob- 
servations with R- and Jf-band photometry from the 
VLT clearly shows that the putative host is quite red; 
R — K = 3.9 mag (V ega). VLT X-shooter observations 
(jKruhler et al.ll2012b[ ) have also recently established the 
redshift of this source based on several strong nebular 
emission lines at z = 2.527. 

The fit to the host SED converges to a young, strongly 
obscured stellar population with a very high current star- 
formation rate of ~ 140M©yr _1 . While rest-frame op- 
tical and NIR observations of this source are limited to 
a single X-band point, its relative faintness (given the 
bright, steep UV continuum) indicates a relatively young 
age and modest stellar mass (~ 4x 10 9 M©). The galaxy 
has characteristics of a luminous, dusty starburst. 

5.2. GRB 051008 

GRB 051008 is the subject of the detailed study of 
Volnova et al. (2013, in prep.), and most of the observa- 
tions reported here will also appear there, although the 
analysis reported in this paper is independent of that 
work. The source was observed within lhr by the 2.6 m 
Shajn telescope at the Crimean Astrophysical Observa- 
tory and by the Tautenburg 1.34 m telescope; although 
early reports suggested some possible optical variation 
of a galaxy pair inside t he original XRT error circle 
(|Rumvantsev et aD 2005a. b), this was not confirmed by 
later analysis, and the most recent XRT error circle ex- 
cludes these galaxies. The large JVh excess rules out a 
high-redshift origin. 

We have imaged the field using LRIS on Keck in numer- 
ous optical filters (U, B, g, R, I, RG$5o) as well as in the 
K band with NIRI. A single, relatively bright galaxy is 
detected consistent with the XRT position in these data. 
The SED of this object shows the unmistakable features 
of a classic Lyman-break galaxy: a flat, steadily rising 
spectrum through most of the optical with a sharp cutoff 
in the B and U bands. Volnova et al. estimate a photo- 
z of 2. 85l°;n5 from these da ta; our own estimate using 
EaZy (jBrammer et alj|2008j) results in z = 2.96 ± 0.11. 
(In principle, this feature could also be fit as originating 
from a Balmer break at z ~ 0.3, but this would require 
a non-star-forming host with very low stellar mass, simi- 
lar to a local dwarf spheroidal. This would be extremely 
surprising for a GRB host galaxy.) We also acquired 
spectra of this source and do not detect any emission 
lines at the location of the host. At the putative red- 
shift only Lyman-a is expected to lie within our spectral 
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range, but because of the resonant scattering properties 
of the Lyman-a line, the nondetection is not a particu- 
larly strong constraint on the p roperties of the ga laxies, 
nor is it unu sual for LBGs (e.g.. IStark et al.ll2010t ) or for 
GRB hosts (jMilvang- Jensen et al.ll2012fl . 

The strong detections in many optical filters allow us 
to constrain the host SFR fairly securely (~ 80 M Q yr _1 
at z = 2.9). The galaxy is young and has a low mass 
(~ 2 x 10 9 M Q ), and (typical of luminous LBGs) it is 
moderately dust obscured (Ay = 0.8 mag). 

5.3. GRB 051022 

GRB 051022 was a High-Energy Transient Explorer 
2 (HETE-2) burst rapidly followed by a range of other 
instruments, including Swift. Early nondetections at op- 
tical wavelengths quickly motivated deeper searches and 
NIR observations, all of which reported only upper lim- 
its on an optical afterglow; comparison of these limits 
to the X-ray flux unambiguously indicates a very large 
dust attenuation column (Ay > 9.3 mag). However, a 
bright host-galaxy candidate was immediately evident 
even in small-telescope images. This GRB and its host 
have been the s ubject of intensive studies by several 
autho rs already ([Castro- Tirado et~aLl l2007bt iRol et all 
l2007f ). which clearly characterize it as a young, high SFR 
ongoing merger at z — 0.809. It seems to have a rela- 
tively high metallicity (Graham et alJl2~009l ). 

Our observations of this source include HST imaging as 
well as a single Keck/DEIMOS /-band image from 2005. 
In addition, the two published papers on this object 
give an extensive library of ground-based optical/NIR 
photometry. Unfortunately, there is a large systematic 
discrepancy between the v alues in two published works : 
most photometric points in ICastro- Tirado et al.1 ()2007bl ) 
are on av erage 0.3 mag br ighter than the corresponding 
values in IRol et all (120071). Since both papers use the 
same calibration ~(lHendenl l2005). and the host is fairly ex- 
tended (~ 1" x 2" in the HST imaging) , we suspect that 
the discrepancy originates from use of different apertures 
under different seeing conditions (neither paper states 
what aperture size is used for the host-galaxy photome- 
try, and the extension is significant even in ground-based 
images) or possibly from the use of extended objects from 
the Henden catalog as calibrators (a significant fraction 
of the objects near the GRB position in the Henden cat- 
alog are themselves significantly-extended galaxies). 

We downloaded the riz GMOS-S images from the 
Gemini Science archive and performed photometry with 
a 2.0" radius aperture using only objects with stel- 
lar PSFs from the Henden catalog, transformed to th e 
SDSS system using the equations of Uester et al.l (|2005f ). 
We also acquired the images from th e Danish 1.54 m 
telescope and UKIRT (published by IRol et al.l 120071) 
and redid the photometry using a similarly large aper- 
ture and the Henden calibration (or 2MASS for the 
NIR bands). These reanalyzed points give much bet- 
ter self-consistency and are used in the model in place 
of the numbers in the litera ture. The photometry from 
ICastro-Tirado et all (j2007bf ) is not used in this analysis. 

Consistent with previous work, we find that the host 
of GRB 051022 is shown to be a moderately young 
(~ 700 Myr), massive (2 x 10 10 M yr" 1 ), rapidly star- 
forming (~ 3OM0yr _1 ) galaxy with moderate dust ob- 
scuration (~ 0.7mag). 



5.4. GRB 060202 

Early, moderately deep nondetections of this 
source first reported by the Faulkes Telescope North 
(jMonfardini et all l2006f) motivated the initial identifi- 
cation of this GRB as a dark burst. The host- galaxy 
candidate was first r eported in twilight Keck imaging by 
ICenko et all (|2006aj) : a faint s ource at this position was 
also seen with Gemini-NIR I ([Schmidt et all I2006D and 
UKIRT (jWang et alj I2006T ). While the co mparison of 
UKIR T and Gemini-NIRI if-band images (jWang et all 
2006) showed no obvious fading, a comparison between 
the AT s -band imaging from the ToO observation and our 
significantly deeper late-time imaging of the field (also 
with NIRI using the K' filter) shows significant fading 
and an astrometric shift of about 0.5" between the 
object centroids, clearly indicating that the early-time 
observations were afterglow-dominated. Since we do 
not have late-time observations in the J or H filters, 
the relative contributions of afterglow and host in 
these bands are ambiguous, but even if the flux in 
these bands was entirely afterglow, the color was much 
redder than a typical afterglow or than comparison to 
contemporaneous X-ray measurements would suggest, 
supporting the notion that this source was reddened by 
host-galaxy dust. 

Our late-time imaging of the field includes only Keck 
g- and i?-band images and the Gemini-NIRI AT'-band im- 
age mentioned ab ove. (We also cautiou sly use the J-band 
photometry from ICenko et al.l l2006al pending indepen- 
dent calibration of the field.) In addition, we acquired 
spectra with LRIS on 2006 Sep. 21 (UT) and detect a 
bright, broad emission line at a wavelength of 6660 A. 
No other lines are evident in the spectrum. We associate 
this line with [O II] at z = 0.78: identification with any 
other common nebular line ([O III], H/3, Ha) would place 
other lines further to the blue where they would almost 
certainly be detected. While the [O II] association would 
put [O III] and H/3 in the detectable spectral range, at 
this redshift both lines land on strong night-sky lines 
that are subject to heavy fringing in the original LRIS-R 
CCD. 

At this redshift, we have only one measurement blue- 
ward of the Balmer break and as a result can only con- 
strain the extinction and SFR crudely, but both are rela- 
tively modest (~ lmag and ~ 6M yr _1 , respectively). 
The faint if -band measurement implies a low mass of 
- 1O 9 M . 

5.5. GRB 060306 

GRB 060306 was rapidly followed by several small tele- 
scopes and no afterglow was detected, although most of 
these limits are not particularly cons training. However, 
the NIC-FPS image ([Chen et al.ll2006fl is quite deep. The 
faint source reported as a possible afterglow candidate in 
this Circular is no longer inside the most recent XRT er- 
ror circle and is most likely a foreground galaxy, but we 
downloaded the reduced image and calculated a point- 
source limit on any afterglow (or host) of K > 18.2 mag 
at the current GRB position relative to 2MASS stan- 
dards. This is several magnitudes below even the mini- 
mum {Pox = 0.5) expectation from the X-ray light curve 
and indicates substantial extinction. 

We acquired late-time imaging of this position in many 
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different optical filters at Keck as well as in the K band 
with NIRI; it was also observed as part of the TOUGH 
project. A source is detected in the XRT error circle in 
all of these observations. In imaging conducted on nights 
with the best seeing, the source resolves into two com- 
ponents: a bright southeastern part and a much fainter 
northwestern source. 

A redshift for t he ho st galaxy of 3.5 was reported by 
iSalvaterra et al.l (|2012t ) based on an NIR emission line 
at 16,800 A, which was claimed to be associated with the 
[O II] doublet. Our reanalysis of the spectrum shows this 
redshift to be incorrect: the apparent separation of the 
line into a doublet is actually due to the velocity offset of 
the two components, an effect that can be clearly seen in 
the 2D spectrum as spatial variation in the line center be- 
tween the two. Extracting the components individually, 
the line is too narrow to be [O II] . 

No other unambiguous lines are detected in the X- 
shooter spectrum, so in principle this line could corre- 
spond to any major emission feature — most likely Ha 
at z = 1.55 or [O III] at z = 2.35, as sugges ted by 
Uakobsson et al.l (|2012| ) and lKriihler et al.l (|2012bl ). Pho- 
tometrically, the lower redshift is favored, as we see no 
evidence of any spectral break due to Lyman-a affecting 
the U band, and a fit to the photometry with EaZy pro- 
vides a consistent z p hot = 1.52+o'Jg (l^ 7 )- Furthermore, 
a weak resolved emission feature is detected in the X- 
shooter spectrum at the expected location of [O II] for 
an assumed redshift of 1.55. Together with the photo- 
metric constraints, we therefore consider the association 
with this GRB and its host system at z = 1.55 reason- 
ably secure. We report photometry of the entire host 
complex (based on a single aperture enclosing both ob- 
jects) in Table El 

The host colors are fairly red with no indication of a 
strong Balmer break, and as with several other galax- 
ies discussed previously we infer a dusty, young, rapidly 
star-forming host, with values (at z = 1.55) of SFR 
« SOOMoyr" 1 , Ay « 2.2 mag, and M « 4 x 10 9 M Q . 

5.6. GRB 060319 

The identification of GRB 060319 as a dark burst is 
based al most entirely o n the early WHT if-band obser- 
vation of iTanvir et al.1 (|2006| ). which reported a if- band 
source inside the XRT error circle. While the coordinates 
given in that Circular are actually not in the most recent 
XRT error circle, reanalysis of the WHT image shows this 
source to in fact be coincident with the most recent XRT 
position, as well as with an underlying object detected in 
late-time Keck imaging with LRIS and NIRC. Compar- 
ison to our NIRC K s -band imaging in particular shows 
obvious fading of this source by almost 2 mag since the 
WHT observation, confirming it as the NIR transient. 
The detected flux is already below the (3ox = 0.5 line 
and the redness of the transient is confirmed by deep 
i?-band nondetections at earlier times. 

We acquired spectra with LRIS on Keck I on 2007 July 
18 (UT). A single, broad line is detected in the red part of 
the spectrum at 8096 A that we associate with the [O II] 
doublet at z = 1.172. No other lines are expected to be 
detected in our spectral range at this redshift. 

A break is evident in the broad-band SED between 
the R and Z bands, consistent with a strong Balmer 



break at this redshift. Actually, the break is sufficiently 
strong that our simplified model of a constant star- 
formation history does not provide a good fit (x 2 /dof 
= 11.9/2), and significant improvement (to x 2 /dof = 
2.7/1) is achieved by allowing a step- function change in 
the SFR to be fit in the model. This causes the fit to 
converge to a current SFR of nominally zero (with the 
UV flux coming from the same older stellar population 
that produces the strong Balmer break) with a 2a upper 
limit of 0.6 M Q yr _1 . While we do not expect the current 
SFR to be actually zero given the detection of [O II] flux, 
the SED is suggestive of a significant decline in the SFR 
since its peak. The mass is well constrained regardless 
of the star- formation history at M = 7 x 10 9 M . 

5.7. GRB 060719 

The position of this GRB was observed rapidly (within 
2 min) by several robotic telescopes, none of which re- 
ported a detection. Deep optical imaging was conducted 
with the VLT (using FORS2) after only 38 min, and 
deep VLT NIR imaging was acquired several hours later 
(using ISAAC). Only the NIR observations resulted in 
the detection of an afterglow. Provisional photometry 
(Malesani et al., private communication) plotted in Fig- 
ure [2] shows that the NIR colors of the transient are 
quite red; our extinction fit indicates a total extinction 
of Ay m 3 mag. 

Photometry of the host galaxy comes from the 
TOUGH project at the VLT and from our Spitzer ob- 
servations. A faint, red object is detected at the after- 
glow position in all of these images. X-shooter spec- 
troscopy iden tified an emission line corresponding to Ha 
at z = 1.532 (Kru hTer et al.ll2~012bf ). consistent with the 
color in the two IRAC channels. The IRAC measure- 
ments and the VLT if -band point constrain the stel- 
lar mass well; we estimate M — 1.3 x 10 10 M Q . The 
SFR and extinction are effectively degenerate parameters 
given the sparse optical photometry (one measurement), 
precluding any independent constraint on the extinction 
from the UV slope and permitting only an upper limit 
on the host extinction to be placed of Ay < 2 mag, but 
for most reasonable estimates of Ay the SFR is relatively 
modest (< 40M Q yr~ 1 ) in this case. 

5.8. GRB 060814 

Bright GRB 060814 was observed at early times by 
a number of telescopes, inclu ding the small ro botic 
TAROT within a few minutes (iKlotz et al.l I2006T) and 
the V LT in the R band after one hour (|Malesani fc Patatl 
2006). No counterpart was reported, although a complex 
extended source near the XRT error circle was noted by 
iMalesanil (120061). Th is source is also detected by SDSS 
(|Ofek fc Cenkof l2006) and no evidence of variability was 
reported. Additional observations several hours later in 
the NIR with UKIRT and P200 detected a K w 18 mag 
source consistent with this location wh ich faded in sub- 
sequent exposures (|Levan et al.ll2006b[ ). confirming it as 
the afterglow of the GRB. 

Our first late-time imaging of this position was car- 
ried out with LRIS during a night of extremely poor see- 
ing (1.8"). The extended source reported in the GCN 
Circulars was still easily detectable (and still obviously 
extended despite the poor seeing); we obtained spectra 
several months later and reported a redshift of 0.84 for 
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this object based on the detection of several nebular lines 
(jThoene et al.ll2007t ). 

We reobserved the source in imaging mode with NIRC 
two years later (in the J and K bands) , and it was also in- 
cluded as part of the TOUGH project (R and K bands). 
The seeing conditions for all of these observations were 
much better (< 0.6"), and these images show that the ex- 
tended object has significant substructure. In particular, 
it is resolved into two distinct subcomponents separated 
by about 1": a point-like object in the southwest, and 
an extended source (consisting of two connected blobs in 
the north and east). 

Additional spectroscopy of this object with FORS-2 
and X-shooter on the VLT was conducted under the 
TOUGH program jjHiorth et all [201 2t Uakobsson et al.l 
l2012t iKriihler et al.ll2012b| ). These observations clearly 
show that the two objects are at different rcdshifts: the 
z = 0.84 emission lines correspond only to the southwest 
object, while the northern/eastern components are at a 
common redshift of 1.923 based on the detection of sev- 
eral emission lines in the NIR. We reanalyzed our LRIS 
spectra and confirm that the z — 0.84 emission lines cor- 
respond only to the western source, with no optical emis- 
sion lines detectable from the eastern source. Relative as- 
trometry between the UKIRT images and the late-time 
VLT images shows that the GRB afterglow is associated 
with the z = 1.92 object, indicating that the z = 0.84 ob- 
ject i s an unrelated foreground system (Pakobsso n et al.l 

The two unrelated sources are heavily blended in our 
Keck/LRIS and Spitzer images. However, given the well- 
determined positions and morphologies of the two objects 
from the VLT imaging and the brightness of both objects 
in all bands, it is not difficult to isolate the flux of the 
true host from the foreground object using the same basic 
technique as for other blended Spitzer images ( §3.5j) using 
a GALFIT model. The photometry reported for this 
object in Table [3] is for the host alone. 

Given the galaxy's remarkable brightness despite its 
redshift of 1.92, it is no surprise our models indicate 
that the host is extremely luminous and rapidly form- 
ing stars. Our best-fit model converges to an SFR of 
~ 240M Q yr~ 1 with an attenuation of Ay — 1.2 mag; 
the mass is moderate (1O 1O M0) and the inferred age 
very young (30Myr), indicating an extremely rapid and 
powerful starburst. This starburst could be coupled to 
ongoing merger activity in light of the apparent binary 
morphology of the host system. 

5.9. GRB 060923A 

GRB 060923A w as the subject of the study by 
iTanvir et al.1 (|2008b| ) , who report deep optical limits and 
a faint K-hanA detection of the afterglow (establishing 
the afterglow as dust obscured) and a basic characteriza- 
tion of the host galaxy. We als o previously discuss ed this 
object as part of the study of iPerlev et al.l ([2009) . Our 
reanalysis here includes all of these data in addition to 
further, unpublished photometry from Gemini-N/NIRI 
and Spitzer (cold mission; PI Fox). 

While the host has been observed in a large number 
of (mostly optical) filters, the detection is marginal in 
all of them. The clearest detections come from the Keck 
£>-band and the VLT i?-band observations; inspection 
of these images hints at a binary morphology with dis- 



tinct components in the northwest (where the afterglow 
is located) and southeast. Given this possible complex 
structure, we have redone the aperture photometry on 
the VLT images to ensure that the same aperture is used 
as for the other optical filters. Our galaxy aperture is 
positioned to attempt to capture most of the apparent 
flux of the possible host complex, but higher signal-to- 
noise ratio (S/N) imaging will be needed to determine 
the exact physical limits of the system. 

Tanvir et al. noted that the host galaxy has an or- 
dinary R — K color, although relative to most GRB 
hosts it is actually fairly red (especially when the Spitzer 
data and measurements in bluer optical filters are con- 
sidered). The redshift of this source is not known spec- 
troscopically: a deep FORS2 integration resulted in de- 
tectio n of continuum down to 4600 A but no emission 
lines (jTanvir et al.ll2012t Ijakobsson et al.ll2012[) . Never- 
theless, the redshift can be constrained by our broad- 
band data; the 95% confidence redshift range from EaZy 
is 1.92 < z < 3.75. The best-fit value is 3.32 with a 
marginalized average redshift of z = 3.05. 

Despite the uncertainty in the redshift, z > 2 is 
strongly favored by the photometry, and the the bright 
Spitzer detections therefore point clearly toward a large 
stellar mass (M ~ 5 x 1O 1O M at z = 2 up to M 
1.5 x 10 11 Mq at z = 3). The derived optical SFR is 
not strongly dependent on the redshift and is high but 
not remarkably so (~ 80 M Q yr" 1 ); given its mass this 
value is in fact quite modest, indicative of an ordinary 
(but massive) high-z galaxy not undergoing a starburst 
episode. 

5.10. GRB 061222A 

GRB 061222A is one of the most obscured afterglows 
in the sample. The only reported detection in the op- 
tical/NIR bands is an early-time Gemini/NIRI if-band 
point, even though this GRB was among the brightest 
events of that year in X-rays and gamma-rays. Deep ob- 
servations in all other filters resulted only in upper limits, 
although these are not as constrain i ng as the -ftf-band 
point. As reported by IPerlev et al.l (|2009h . the host is 
quite blue, and in fact shows a strong Lyman-a emission 
feature. It i s the only dark GRB host k nown to exhibit 
this feature (jMilvang- Jensen et al.l[2012T ). 

In addition to a reanalysis of the photometry presented 
in that work, here we include several new observations: 
specifically, Spitzer/IRAC observations and an archival 
FST/NICMOS F160W measurement (PI: Berger). Con- 
sis tent with the very blue color of this object reported 
by IPerlev et al.l (|2009l ). the Spitzer data only marginally 
detect the host (it is one of only a few sources in the 
sample which do not have a strong Spitzer detection), 
but the source is well detected by HST. 

Given the very blue rest-frame UV color (which may 
indicate an even younger starburst than our 5 x 10 7 yr 
minimum value) the extinction is extremely low, converg- 
ing to the physical limit of Ay — for nearly all trials, 
even if a significantly lower metallicity is assumed (the 
extremely blue color and low mean extinction mark the 
host of GRB 06 1222 A as a dramatic outlier compared 
to most other dark GRB hosts in our sample). This is 
particularly curious given the tremendous extinction evi- 
dent in the afterglow SED, and suggests extreme hetero- 
geneity in the galaxy's internal dust distribution: either 
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the dust is localized to a small number of very optically 
thick clouds and the GRB sightline happened to inter- 
sect one, or the bulk of the star formation is actually 
heavily obscured. The small stellar mass of this object 
probably favors the former interpretation (the estimated 
stellar mass of ~ 10 9 M© is much lower than that of typ- 
ical submillimeter galaxies), but deep long-wavelength 
observations (e.g., with ALMA or EVLA) will be needed 
to resolve the question unambiguously. 

Interestingly, the nearby galaxy seen in our images at 
an offset whose relative orientation and proximity might 
suggest a merging companion is not at the same redshift 
(|Perlev et al.ll2009f ). 

5.11. GRB 070306 

GRB 070306 was the subject of Uaunsen et all (|2008t ). 
who identify it as a highly obscured afterglow based on 
deep, early-time WHT and VLT follow-up observations: 
the afterglow is detected only in the K and H bands 
in these data, and it is very red. Data points for our 
modeling (and the host redshift of 1.496) are derived 
from this work (Tab le 2) and from the recent study of 
IKruhler et al.l (|201lD . supplemented by our Spitzer pho- 
tometry. This combined dataset is one of the most com- 
plete of any object in the sample. 

The host galaxy is relatively bright and well detected 
in most filters and has a relatively blue SED showing a 
clear Balmer break that matches the emission-line red- 
shift. The stellar mass and SFR are both substantial 
but the mean extinction is relatively low. The fit quality 
is actually improved significantly by allowing a noncon- 
stant star-formation history that is elevated in the past 
relative to the present; even with this flexibility added, 
our model converges to a well-determined SFR (a rela- 
tively modest ~ 7MQyr _1 ) and significant stellar mass 
(~4x 10 10 M ). Given the blue rest-frame UV slope the 
extinction is very low, Ay < 0.3 mag. 

As with GRB 06 1 222 A (and the recent GRB 100621A; 
IKruhler et alJ[20TlT l. the combination of blue host and 
highly reddened afterglow implies either a very hetero- 
geneous ISM or a highly obscured super-starburst occur- 
ring inside the galaxy core. The much larger mass of this 
galaxy is consistent with either picture. In fact, we have 
recently detected this object at 2 GHz with the EVLA 
(Perley et al. 2013b, in prep.), suggesting that in fact 
this object is indeed a blue submillimeter galaxy (the in- 
ferred SFR from the radio continuum is ~ 200 M yr _1 .) 

5.12. GRB 070521 

GRB 070521 is an extremely dark burst; it occurred 
during dark time and triggered both the Faulkes Tele- 
scope North and P60 for rapid observations, neither 
of which detecte d a transient (|Melandri et al.l l2007t 
ICenko et al.ll2009T ). Follow-up observations within 24 hr 
were conducted with P200, Gemini, Keck, and Subaru, 
none of which reported any variation. (Some possible, 
ve ry marginal evidenc e of optical variation was suggested 
bv ICenko et~aT1 120071 and IXin et al.ll2007l but neither of 
the reported positions are in or near the final XRT error 
circle, nor would variation in these optical bands be ex- 
pected given the much deeper NIR limits available from 
Gemini.) Given the deep, early NIR limit and bright X- 
ray afterglow, the inferred limit on the extinction column 
is extremely large (Ay > 12 mag). 



The best position available for this source is the 
UVOT-enhanced XRT error circle. There is only one 
source that is clearly inside this circle: a very red, ex- 
tended object toward its eastern side. A fainter (in most 
bands), bluer, morphologically complex source is located 
immediately southeast of this galaxy, just outside the 
XRT error circle. It is not completely unambiguous 
which of the two sources is the true host galaxy (or if 
they are related), but given the first object's brightness 
and closer proximity to the error circle we identify it as 
the host galaxy. This is the same object w e previously 
identified as the host in lPerlev et al.l 120091 (our imaging 
was not deep enough to identify the fainter galaxy in that 
work). 

We previously estimated (|Perlev et al.l 12009) a photo- 
metric redshift of 1.35 for this host. Our additional ri- 
band images and a spectrum of the host which identified 
no sign of [O II] emission out to about 10,000 A now lead 
us to favor a slightly higher redshift of about 1.7 (the 
95% confidence range from the EaZy fit is 1.37-2.20). 
We continue to see evidence of a strong Balmer break 
blueward of the J band that suggests an evolved stel- 
lar population and significant mass — in this case about 
3 x 10 10 M Q . The SFR is modest at ~ 40 M Q yr" 1 , and 
the galaxy is also quite dust obscured {Ay ~ 2.2 mag). 

5.13. GRB 070802 

The relatively faint, red afterglow of this event was 
observed by several telescopes including GROND on 
the MPI/ESO 2.2 m, which detec ted the afterglow i n 
all seven filters from g though K (jKriihler et al.l 12008). 
(GRBs 080607 and 081109 are the only other events in 
this sample to be detected in the R band or bluer fil- 
ters.) The optical afterglow was also observed spectro- 
scopically wit h the VLT, establishin g an absorption red- 
shift of 2.454 (jEhasdottir et alJl2009T ). Both the GROND 
SED and the VLT spectrum show a clearly reddened 
continuum with a significant detection of the 2175 A 
dust bump. The actual degree of dus t atten uation is 
borderline for inclusion: IKruhler et al.l (|2011h estimate 
Ay = 1.23 ± 0.17mag, only marginally ruling out the 
Ay = 1 mag threshold. 

Most of our photome tric data points are tak en from 
the VLT observations of Ehasdo ttir et al.l (I2009T) and the 
GROND observations of IKruhler et all (|2008t l. We add 
to this our Spitzer observations of the source from our 
program. The host position in the 3.6 /xm image was 
affected by column pull-down originating from a bright 
star elsewhere in the image; we removed this effect by 
adding flux to this column in the vicinity of the host until 
the artifact disappeared and then performed photometry 
in the usual way, but this point should nevertheless be 
used with caution. 

While well detected in the optical bands, the host of 
GRB 070802 is only weakly detected by Spitzer, indicat- 
ing a low stellar mass and relatively young population 
(although the sparse optical data make the actual SFR 
and dust attenuation difficult to characterize). We esti- 
mate a SFR of ~ 8OM yr _1 (uncertain to a factor of 
~ 2) and a mass of ~ 10 9 M®. 

5.14. GRB 071021 

GRB 071021 was a very long burst initially flagged as 
a high-redshift candidate based on its prompt emission 
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properties (|Sakamoto et al.l I2007D . a possibility seem- 
ingly encouraged by the lack of detectio n of an op- 
tical counterpart in early imaging (e.g., iHsiao etaLl 
120071: iGuidorzi et al.l |2007| ). Subsequent imaging did 
detect a faint, fading NIR counterpart with red 
KJz colors consistent with a dust-reddened afterglow 
(jCastro-Tirado et alj|2007al) . 

We have imaged this field in several different filters 
with LPJS and NIRC on Keck, as well as with NIRI (J 
band during cloudy conditions). A faint galaxy is de- 
tected at the afterglow position in all of these images 
except J. The surrounding environment is complex; the 
host galaxy is part of a "chain" of four different sources 
along a line separated by a few arcseconds. The associ- 
ation of these galaxies with each other is not clear (the 
nearest source, about 2" southeast of the host, has a 
similar, red color), but fortunately the host is well sep- 
arated from the other objects and is relatively bright, 
allowing construction of a good SED. This source was 
al so included in the X- shooter spectroscopic campaign 
of iKruhler et al.l ()2012bl ) and has a measured redshift of 
2.452. 

The combination of bright Spitzer fluxes and bright 
optical detections indicate a massive and rapidly star- 
forming galaxy. The SFR is very high (SFR w 
100 Mq yr _1 ), as are the stellar mass (~ 10 10 M©) and 
the extinction (Ay ks 1.6 mag). The J nondetection just 
blueward of the Balmer break is constraining and sug- 
gests a significant population of evolved stars, although 
the host is still fit well by a continuous star-formation 
model. 

5.15. GRB 080207 

GRB 080207 was followed quickly with deep imaging at 
the VLT and not detected in any filter, including in the 
NIR. The corresponding upper limit on the extinction 
(Ay > 2.0 mag) is actually not as strong as that of most 
other sources in our sample since the X-ray afterglow 
faded quite rapidly, but the depth of the nondetections 
(coupled with a very high X-ray Nn) clearly identify this 
event as a dust-extinguished dark burst. 

The host galaxy of GRB 080207 has been studi ed in 
detail in previous works by ISvensson et al.l (|2012T) and 
iHunt et all (|2011h . The photometric redshifts presented 
in these two papers (z — 1.74+ ;jg and z — 2. 2+ ;g, 
respectively) have also recently been confirmed by X- 
shoote r spectroscopy that establishes the redshift of 
2.086 (Kriihler et al. 2012b). Our analysis here combines 
photometry from these two sources and fixes the redshift 
to the spectroscopic redshift. (Except in the case of the 
NIRC i^-band observation, we have not reanalyzed the 
photometry independently.) 

Because the source is so red (in terms of R — K it is the 
reddest object in the sample), and the optical detections 
are all at relatively low significance, there is significant 
degeneracy between a mature population and a highly 
dust-obscured, rapidly star-forming one. This results in 
a large range of uncertainty in the SFR and Ay , and even 
an SFR of zero is only ruled out at about la, although the 
clear detection of several emission lines in the NIR with 
X-shooter suggests that the actual value is probably not 
far from the best fit, ~ 50 M Q yr _1 (the dust- uncorrected 
Ha SFR from the X-shooter provides a lower limit of 



10-15 M Q yr" 1 ). In addition, the HST magnitudes are 
puzzlingly discrepant (by about ~ 0.5 mag, many times 
the photometric error) from ground-based photometry 
in the nearby filters, leading to an unavoidably poor fit. 
The mass is not affected by these issues and is well con- 
strainted at ~ 1.2 x 10 11 M©, a large value intermediate 
between (but co nsistent wi th) the estimates p rovided by 
ISvensson et~aT] (|2012f ) and lHunt etall (|2011f ). 

5.16. GRB 080325 

This e vent and its host galaxy were the subject of the 
study of lHashimoto et alT(|2010l ). A faint, red, fading 
NIR transient was detected by Subaru at a position off- 
set slightly from a fainter, extended source that repre- 
sents the presumptive host of the GRB. The photomet- 
ric redshift of the galaxy from that work is z = 1.9^ i- 
Further Subaru observations since that time using the 
MOIRCS spectrograph have confirmed the redshift of 
this host to be 1.78 (Hashimoto et al., private commu- 
nication). Our own EaZy photometric redshift estimate 

(%>hot — l-57io 15) i s a ^ so consistent with this value. 

We have observed this field with LRIS in several opti- 
cal filters and with Spitzer. These observations confirm 
the characterization of this source as a faint, red galaxy 
(although we measure s omewhat different optic al magni- 
tudes than reported by lHashimoto et al.ll2010h . 

Given the very bright K-bemd and Spitzer fluxes 
this galaxy is quite massive. Our step-function star- 
formation history initially converged to a maximum stel- 
lar age older than the Universe at this redshift; to prevent 
this unrealistic behavior for this event we modified the 
model to an exponentially falling star-formation history, 
which produces a cosmologically consistent answer. The 
SFR of this host is relatively modest (~ 15 M Q yr _1 ) in 
comparison to its total mass (~ 1.5 x 10 10 M Q ). These 
values are broadly con sistent with those reported by 
lHashimoto et al.1 (20101 ). 

5.17. GRB 080607 

This exceptional event has been the subject of sev- 
eral previous papers by our collaboration detailing vari- 
ous aspects of the burst and its host environment. No- 
table attributes of this burst include the e xtreme lumi- 
nosit y of its prompt emission and afterglow (jPerlev et al.l 
|2011| ). well-determined line-of-sight extinction proper- 
ties (Ay « 3 mag of dust wit h a clear 2175 A absorp- 
tion bump: iPerlev et al.l 1201 1^ . an afterglow spectrum 
rich in ionic and molecular absorption features from a 
molecular cloud within the host along th e line of sight 
(jProchaska et al.l 120091 : ISheffer et aHl2009f ). and the rel- 
atively high redshift of the system (z = 3.038 from 
the same spectra). The ho st galaxy was previously 
studied by [Chen et al.l (|2010f ) ( see also the Erratum at 
iChen et al.N2011h . 

We have reanalyzed all of our ground-based photome- 
try of this event for the study here. By fixing the aper- 
ture at the host location we achieve marginal (2-3cr) de- 
tections of the host in both the Magellan r an d Keck / 
bands. The basic result of IChen et aH (|2011h remains 
the same, which is that this is a massive, star-forming 
host. Compared to other galaxies in our sample, how- 
ever, none of its properties are particularly extreme (SFR 
wl8M yr 1 ,Affti4xlO lo M Q ). 
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5.18. GRB 081109 

GRB 081109 is an event analogous to GRB 080607, 
a luminous burst occurring within a dusty environment 
which heavily reddened (but did not prevent detection 
of) its optical afterglow. No spectra of the afterglow were 
available, but the redshift 0.979 was established from 
later spectroscopy of its bright host galaxy (the location 
of which is coincident with the afterglow based on com- 
parison of early-time and late-time GROND imaging). 
While none of the optical measurements are below the 
/?ox < 0.5 criterion, the GROND afterglow SED provides 
a strong constraint on the extinction of Ay = 3.4+q 3 
(|Kriihler et all |2~011| ). easily satisfying the criterion for 
inclusion in this sample. 

Our observations are taken mostly from Kriihlcr ct al. 
(| 2 llf) and supplemented by our Spitzer photometry, 
which provides an improved (slightly higher) estimate 
of the stellar mass. Our va lues are otherwise gener- 
ally consistent with those of iKruhler et al.l 120111 (SFR 
= SOMgyr" 1 , M = 9 x 1O 9 M , Ay = 1.3 mag) and 
indicate a moderate-mass, rapidly star-forming host. 

5.19. GRB 081221 

Another very bright Swift burst, GRB 081221 was not 
observed rapidly, but deep observations with a variety 
of telescopes at later t i mes reported no op tical afterglow 
(|Malesani et al.l 120081: lAfonso et all 120081) with the ex- 
ception of NIR imaging with Gemini-Nort h, which de- 
tected a faint source in only the K band (|Tanvir et al.l 
l2008aT) . This source was observed to fade in our later 
Gemini epochs, leaving a slightly extended source un- 
derlying this position. The source is also visible in opti- 
cal images (where it do es not fade) , a s also reported by 
iPerlev et all pOOl and lAfonso eTall pOOl . The region 
around the transient is fairly complicated, with several 
other galaxies visible within a few arcseconds. None of 
these appear to be visibly interacting with the GRB host. 

Photometry is reported in a variety of optical and NIR 
filters from Keck and Gemini observations. Our late- 
time NIRI imaging in the J and H filters occurred during 
cloudy conditions and are less constraining than taking 
the 2a upper error bar on the early-time measurements 24 
as an upper limit on the host flux, so we use the latter 
in our modeling and in Figure [3] This galaxy was also 
observed with X-shooter; a secure redshift based on de- 
tection of the [O III] and Ha emi s sion l ines at z = 2.260 
was reported by ISalvaterra et all (|2012l ). 

The host of this system is massive (~ 6 x 10 10 M Q ) but 
young and rapidly star forming (~ 70 M Q yr" 1 ): our H- 
band limit in particular suggests a fairly subdued Balmcr 
break given the well-determined optical flux and slope. 
As with all other rapidly star-forming hosts in our sam- 
ple, the dust extinction is also large (Ay — 1.4 mag). 

5.20. GRB 090404 

GRB 090404 was imaged rapidly by the Xinglong TNT 
and after about 8hr with the NOT and (in the NIR) 
the Calar Alto telescope; no transient was detected with 
any of these observations (|Xin et al. 2009; Malcsan i et all 

24 These likely contain some unquantified afterglow contribu- 
tion; unfortunately, late-time observations are not yet available for 
subtraction. 



120091 ; Uelanek et all l200"9l) . Howev er, an afterglow was 
detec ted at millimeter wavelengths (jCastro-Tirado et~all 
f2009h . 

Our imaging includes Keck B, g, and / bands, as well 
as HST and Gemini/NIRI {K band) taken 22 days after 
the burst. The position of the afte rglow is established 
by th e millimeter observations of ICastro-Tirado et all 
(2009), which place it near the center of an extended, 
highly elongated source that is clearly detected with HST 
(in F160W only) and NIRI. It is marginally detected in 
most optical filters. The field around this presumptive 
host is complex: a pair of bright interacting galaxies is 
evident slightly to the east (the near edge of one galaxy's 
disk is only 1.5" away) and a compact, blue source is 
visible 1" to the northeast. These nearby sources po- 
tentially can introduce significant blending; fortunately, 
all images reported here were taken in very good see- 
ing conditions, and the most nearby source is very blue 
and contributes negligibly to the Spitzer bands. (Never- 
theless, as with all other Spitzer observations in crowded 
fields, we subtracted all nearby sources before performing 
photometry) 

No spectroscopic redshift is currently available for this 
object. The galaxy is consistent with a photometric red- 
shift of ~ 3 based on a (questionable) Lyman break in 
our photometry; we do not significantly detect the host 
in B despite a relatively strong detection in g. Although 
the -B-band image is slightly shallower than the g-band 
image, this nondetection is about 1 mag below the ex- 
pectation from an SED model extrapolated from the 
other filters. This redshift is generally consistent with 
the other colors available for the object, in particular 
the suggestion of a Balmer break from the I-H-Z color 
and the relatively red IRAC 3.6-4.5 /zm color. However, 
this association is not highly significant (EaZy gives only 
1.25 < z < 3.84 at 95% confidence, or z = 2.87±g'ls at 
la). 

Assuming the Lyman-a break is real and the red- 
shift is high, the bright IRAC detections indicate a 
quite massive galaxy: ~ 10 11 M at z = 3, simi- 
lar to the values measured for GRBs 080207, 080325, 
and 060923A. The extinction is significant but relatively 
modest (Ay = 1.1 mag) and the SFR, while not well 
constrained as a result of the low-S/N optical detection 
rates, is ~ 40 M Q yr _1 ). 

5.21. GRB 090407 

Early-time observations of this GRB were limited to 
the small ROTSE telescopes (which did not detect a 
cou nterpart but are not d eep enough to be constrain- 
ing; ISchaefer fc Yua n 2009) and the X-ray afterglow was 
relatively faint at this t ime, but deep VLT i?-ba nd and 
z-band limits at 13.7 hr dMalesani fc Fvnboll2009l ) are far 
below the extrapolated X-ray flux and identify this event 
as a da rk burst. GROND limits (including NIR obser- 
vations; Afonso et al. 2009) were contemporaneous with 
the VLT observation but, given the time since the burst, 
these are not particularly constraining. 

The host galaxy of this event is extremely red, and to 
date we have clearly detected it only in the NIR, with 
HST, Spitzer, and Gemini-NIRI (despite cloudy condi- 
tions) , although forced aperture photometry at the after- 
glow location does give weak (3<r) detections in our Keck 
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obser vations. An X-shoo ter redshift of 1.448 is avail- 
able (|Kriihler et alJl2012b[ ) based on detection of Ha and 
[O II]. Given the weak (and possibly uncertain) optical 
detections the properties of the system are difficult to 
constrain with any precision (our parameter estimates 
are all uncertain to an order of magnitude), although 
as with most other galaxies in our sample, moderate- 
to-large values for the SFR, extinction, and mass are 
preferred. 

5.22. GRB 0904 1 7B 

GRB 090417B was a moderately bright but very 
long burst; its (bright) X-ray afterglow was immedi- 
ately recognized to be coincident with an SDSS galaxy 
(|Fvnbo et al.ll20p9af) whose redshift was measured early 
on to be 0.345 (jBerger fc Foxl [20091 ). Given this bright 
host, afterglow upper limits reported in the GCN Circu- 
lars are not straightforward to interpret (image subtrac- 
tion against a later reference epoch would be necessary); 
in the very conservative assumption that the afterglow is 
no brighter than the quoted host-galaxy magnitudes, the 
Ay > 1 mag criterion is marginally satisfied, though in 
all likelihood the real limits are much deeper than any- 
thing currently quoted and so the actual extinction is 
probably significantly larger. No thorough study of the 
optical afterglow has yet been performe d; the claimed 
WOT limit of u > 24.9 mag at 11 hr in IHolland et al.1 
( 2010) does not seem reasonable given typical limits at- 
tained by the UVOT in a few hours of exposure time. 
(The specific estimate of Ay > 11 mag in IHolland et al.1 
120101 is based on the Nn measurement from the X-ray 
afterglow, not on optical limits.) 

As the lowest-redshift event in the sample, the data on 
the host galaxy of this extremely obscured event are par- 
ticularly extensive, including several UVOT detections 
(not shown in Figure [3]) as well as a wide suite of detec- 
tions through the optical-NIR bands and Spitzer pho- 
tometry. The major ity of our measurements come from 
IHolland et al.l (|2010l) , and are supplemented by our HST 
and Spitzer photometry 

Like some other bursts in the sample, the strong 
Balmer break points toward a declining SFR despite a 
relatively young stellar age. For a constant SFR we de- 
rive a best-fit value of 3 M Q yr _1 , but a significantly bet- 
ter fit is obtained if the current SFR is decoupled from 
the assumption of a constant star-formation history (in 
which case it converges to 0.35 ± 0.2 M Q yr _1 , suggesting 
a recent decline in the SFR). The stellar mass is moder- 
ate at ~ 3.5 x 1O 9 M0, intermediate between the SMC 
and LMC. Interestingly, while the host of GRB 090417B 
is a relatively low-mass galaxy (like other low-z GRB 
hosts), it does not appear to be forming stars at a par- 
ticularly high rate for its mass (quite unlike most other 
low- z GRB hosts). The host is dusty {Ay = 0.9 mag) ac- 
cording to our fits, although this value is lower than for 
other events in our sample and is far less than the value 
of Ay = 3.5±i° quoted by IHolland eTaLll20lol which as- 
sumed an intrinsic power-law host SED (an assumption 
which is unreasonable for this galaxy given the strong 
Balmer break). 

5.23. GRB 090709A 

GRB 090709A was another extremely bright Swift 
burst and had an extremely bright X-ray afterglow. The 



event triggered both P60 and PAIRITEL, which auto- 
matically slewed to the location and detected only a very 
faint and red afterglow. Our analysis o f the af terglow was 
previously published by ICenko et al.l (|2010l ). in which 
we interpret the burst as a classical dark burst with a 
highly dust-reddened afterglow. However, the possibly 
quasiperiodic prompt-emission behavior of this event led 
some to suggest that the event may actually be a Galac- 
tic source such as a magnet ar (Ma rkwardt et al.l 120091 : 
iGotz et al]|2009f) . This claim has been disputed by oth- 
ers on account of its high Galactic latitude, evidence for 
X-ray absorption, and the classical behavior of its after- 
glow (|de Luca et al.ll2010t ICenko et alj|2010t ). 

Host-galaxy observations reported in our earlier work 
included deep integrations of the field with GTC and 
Keck (LGS-AO) which did not detect a host galaxy at the 
afterglow position. Since that time we have integrated 
much deeper in both the optical (with LRIS) and NIR 
bands (with Gemini), as well as with HST (WFC3/IR) 
and Spitzer. A faint object is detected in all of these ob- 
servations (except in the i?G850 filter) directly at the po- 
sition of the NIR transient. In the HST observations this 
object is clearly elongated in the E-W direction, identi- 
fying it as a distant galaxy. 

The HST photometry, combined with the deep i?G850- 
band limit from LRIS, clearly indicates the presence of 
a spectral break in the SED at around 11,000 A which 
we associate with Balmer absorption at z « 1.8. A fit to 
the broad-band photometry with EaZy gives a 95% con- 
fidence redshift ra nge of 1.14 < z < 2.3 4, which also rules 
out the claim by Ide Luca et al.l (|2010f ) of a moderately 
high redshift of 4.2 ± 0.2 based on the X-ray absorption 
properties. (Indeed, the 5-band detection of the host 
alone rules out this redshift.) 

The properties of the host are relatively moderate 
among events in this sample: at z — 1.7 we have 
M w 10 10 M Q , SFR » SMeyr" 1 , and Ay w 1.4 mag. 
The SFR is quite sensitive to the assumed redshift; if 
the redshift is at the bottom end of the allowed range 
then only the g and B filters measure the rest-frame UV 
(which are too similar to usefully constrain the spectral 
index and therefore the extinction correction to the SFR) 
and the SFR cannot be usefully constrained. However, 
as long as z > 1.3 (which is secure at ~ 90% confidence) 
this is not a concern, and the properties are similar to 
those outlined above. 

6. SAMPLE PROPERTIES 

The properties of the host sample discussed above 
differ starkly from those of the optically selected host 
samples of the past. While most such previous efforts 
have struggled to ev en detect ~ 50% of hosts in the 
rest-frame NIR (e.g.. iLe Floc'h et al.ll200l l2006ft . every 
host in our sample was detected in the NIR and with 
Spitzer (Figures [SHHl and EH5]) . The host galaxies in our 
sample show a diverse range of colors but include two 
extremely red objects (EROs, with R — K > 5 mag: 
GRBs 080207 and 090407) and many additional objects 
that are just short of this threshold (R — K > 4 mag: 
GRBs 090404, 090709A, 080607, 070521, 071021, 060306, 
060923 A, and 080325). Only two GRB host galaxies 
with R — K > 4 mag were known pie- Swift (|Levan et al.l 
l2006al iBerger et al.ll2007t ). both of which were also dark 
bursts (Figure [9]). 
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To facilitate direct comparison between our results and 
optically selected GRB host studies, we downloaded the 
photometric data for all pre- Swift GRBs with known red- 
shift and deep, previously published observations in both 
the optical and NIR bandpasses f rom the online com- 
pilation in the GHosts database 25 (Sav aglio et al.ll2009t ) 
and ran these SEDs through the same procedure as for 
the galaxies in our obscured-burst sample, fitting for the 
SFR, mass, and Ay using our SED analysis code. The 
sample (and the results of our SED fitting) is summa- 
rized in Table [4] Most of these events are optically bright 
bursts, although a few are "dark" and probably similar in 
nature to the objects in our sample. Our derived values 
ar e generally cons i stent with those previously published 
by ISavaglio et al.l (|2009f ) , although we were also able to 
solve for (or at least constrain) Ay for all hosts, which 
Savaglio et al. attempted only for a few objects. We did 
not use any long- wavelength > 10 /im observations: an 
analysis of submillimeter/radio SFRs will be reserved for 
future work. 

We did not extend this comparison sample to Swift- 
era, bursts. Only a few Swift hosts have been studied 
at the level of detail necessary to perform SED fitting 
(3+ broad-band filters, including NIR observations), and 
these events may not be representative since the pursuit 
of multi-filter observations and publication of the result- 
ing work was probably influenced by the properties of 
the hosts in the first place. These same concerns proba- 
bly apply to some extent to pre- Swift hosts as well, but 
since the sample of known-redshift events was relatively 
small before 2005, observations were necessarily concen- 
trated on a limited number of objects (specifically, our 
sample of 31 "well-observed" hosts constitutes more than 
65% of all pre- Swra/f events with known redshifts) — and 
so the literature sample of published host photometry is 
much less susceptible to selective-observing and selective- 
publication bias than in the current era. 26 

Due to the lower sensitivity of most other GRB satel- 
lites relative to Swift, the median redshift of pre- Swift 
events is much lower (z ~ 1—1 . 5) th an for Swift events 
(z « 2-3: Uakobsson et al.ll2006ll2012[ ). As a result, while 
the pre- Swift sample is appropriate for comparisons at 
z « 1, it is much more limited at higher redshifts. While 
there does not yet exist a survey of Swift GRBs with 
both understandable selection effects and deep observa- 
tions in enough filters to model the SED in any detail, 
some individual studies using only one or two filters have 
recently been published that are useful for extending our 
comparisons to higher redshifts in a more limited sense. 
The Optically Un biased GRB Host Project (TOUGH; 
iHiorth et al]l20T2 ) is a large ground-based VLT project 
centered on deep R and K imaging of 69 uniformly se- 
lected bursts that is 80% redshift-complete. (Since this 
sample is not optically selected it includes dust-obscured 
events; indeed, six dust-obscured events in this paper 

25 grbhosts.org . 

26 Still, even the pre-Swift sample is probably not perfectly rep- 
resentative of the hosts of optically bright bursts; a bias against 
very faint hosts may still be present (since acquiring multiple fil- 
ters may not be attempted if the first observation produces a deep 
nondetection, and emission-line redshifts are much easier to pro- 
vide for brighter hosts). As we continue to find an overabundance 
of faint galaxies after incorporating the obscured population, we do 
not expect that this significantly impacts the primary conclusions 
of this work. 
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Fig. 7. — A"-band apparent magnitudes of GRB host galaxies 
from our sample (Ay > 1 mag Swift burs ts, squares) and others: 
pie-Swift GRBs from Savaglio et al. 2009, diamond s; a un iformly 
selected sample of Swift GRBs from Hjorth ct al. (2012), smal l 
circles; and optically bright z > 3 GRBs from Laska r et al.ll20lH . 
large circles. Colored symbols indicate GRB hosts, with dark red 
indicating dust-obscured events. Hollow symbols indicate it-band 
faint hosts (K > 23 Vega mag). Additionally, i n light gray we 
plot fi eld galaxies from the MOlRCS Deep Survey I IKaiisawa et al.l 
|20TO|) of the GOODS-N field, which is complete to K = 23.1 at 5<r. 
The area of each field-galaxy symbol is scaled by the galaxy's SFR, 
and so above the completion limit of the survey (dashed horizontal 
line) the GRB host distribution should follow the distribution of 
gray "ink." We have quantified this by plotting dotted and solid 
curves representing the lst/3rd quartile boundaries and the median 
(respectively) of the SFR-weighted distribution of K < 23 mag 
(Vega) field galaxies and, therefore, the expected distribution for 
K < 23 mag GRB hosts if they follow star formation uniformly. 
This expectation is certainly not met at 2 fS 1; a large majority 
of hosts occur in galaxies fainter than the expected median (but 
brighter than the MODS completeness limit). NOTE: Figures 7- 
13 show only a subset of the MODS galaxies due to arXiv file size 
limitations. 

overlap with the TOUGH sample.) The K-h&nd limits 
in TOUGH only reach K « 21 mag, but a much deeper 
survey of the hosts of a sample of randomly selected, op- 
tically bright Swift bursts at at z < 3 was conducted 
with HSTm the F160W filter (Tibbets-Harlow et al., in 
prep.). Looking to even higher redshifts, a recent IRAC 
study of optically selected GRBs at z > 3 (|Laskar et al.l 
1201 If ) provides deep observations of a large sample of 
high-z hosts at 3.6 /im. 

In addition to comparing the GRB populations to each 
other, a major motivation of this work is to compare 
both populations to field galaxies in order to reevalu- 
ate to what extent the inclusion of dust-obscured GRB 
hosts in the overall sample is likely to alter our per- 
ceptions of how well the GRB rate tracks the SFR at 
high redshift. While a rich literature of deep compari- 
son surveys of many different fields now exists, we have 
chosen as our primary comparison sample the catalogs 
provided from the M OIRCS Deep Survey (MODS) of 
iKaiisawa et al.l (|2010l[20TH ). a small- area, ultradeep sur- 
vey of the GOODS-N field, selected in the fT-band and 
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Fig. 8. — F160W (i/-band) apparent magnitudes of GRB host 
galaxies from our sample and pre- Swift comparison sample, plus 
a sample of randomly selected Swift GRB hosts with measured 
optical redshifts (representative of optically bright Swift bursts) 
at z < 3 from the HST Snapshot survey (Tibbets-Harlow et al., 
in prep.) Symbol conventions are the same as in Figure [7] except 
that we plot the Snapshot hosts (cyan triangles) in place of the 
TOUGH galaxies. These hosts confirm the trend seen in Figure [7] 
at z £S 1 for GRBs to prefer subluminous galaxies, but at z « 2 
both optically selected and dark GRBs commonly occur in lumi- 
nous galaxies and show no strong deviations from the expected 
magnitude distribution. 

complete to K — 23 Vega mag with redshifts measured 
by a combination of spectroscopic and photometric tech- 
niques. The primary driver of this choice of comparison 
is the depth of the MODS survey (reaching well below 
L*) and its selection in a single, mostly dust-unbiased 
NIR band, avoiding the need to extrapolate or make 
large completeness corrections in order to quantify the 
importance of the faint galaxies which many GRBs are 
observed to inhabit. Magnitudes and derived properties 
are taken directly from the online catalogs, with an ad- 
justment to the SFRs and stellar masses of a factor of 1.8 
(decrease) to convert from a Salpeter to a Ch abrier IMF 
in order to fairly compare these results to ours ; Ki'h et ahl 
I2006t iMichalowski et al]|2012al) . 

6.1. Redshift Distribution 

While the redshift distribution of our sample is not 
expected to necessarily reflect the parent Swift distri- 
bution given our selection methods (Sj2|), we actually 
observe the two to be reasonably similar (a K-S test 
against a catalog of all publicly available Swift red- 
shifts gives a probability value of 0.36 that the distri- 
butions are drawn from the same population). This 
suggests that any redshift dependence in the intrinsic 
dust obscuration of GRBs is relatively minor. There 
is a possible hint toward a dearth of relatively nearby 
GRBs: half of the 5wi/£-triggered events (11 out of 
22) are in a narrow redshift range of 1.5-2.5, an epoch 
centered on the peak of cosmic star-formation inten- 



FlG. 9. — R — K color of GRB hosts and field galaxies. Sym- 
bol conventions are the same as in Figure [7J and the curves in- 
dicate quartile boundaries of the expected host distribution for 
K < 23 mag galaxies. The hosts of z > 1 dust-obscured GRBs 
are much redder than those of unobscured hosts. GRBs clearly 
appear to avoid red hosts at 2 a 1 compared to expectation from 
the relative SFRs of this population, although this effect is signifi- 
cantly less pronounced (or may even be absent entirely) at higher 
redshifts. 

sity a nd of the abundance of submillimeter galaxies 
(e.g.. | Chapman et al.1 120051: iHopkins fc Beacoml 120061 : 
iWardlow et al.l 120111: IMichalowski et al.l l2012bl) . Only 
about 25% of all Swift GRBs are in this same range (e.g., 
Uakobsson et al.l 120121) . More observations of a cleaner 
sample would be needed to confirm whether this repre- 
sents a real effect (rather than a more subtle result of 
selection bias or a result of small- number statistics) . 

6.2. Stellar Mass 

The low observed stellar masses of the population of 
known GRB host galaxies has long been one of the 
strongest arguments for a GRB population that does not 
trace star formation uniformly in all environments. Many 
z w hosts have extremely low mass, falling well below 
the stellar masses of typical star- forming galaxies selected 
from SDSS or from the hosts of Type II or Type Ib/c 
super novae found in untargeted surveys (|Modiaz et al.l 
2008). While it is not completely clear how the low- 
redshift (z < 0.1), underenergetic (Ei so < 10 49 erg) GRB 
population that dominates these samples relates to more 
typical 27 cosmological GRBs, generally (if not ubiqui- 
tously) low masses are also observed in the case of the 
host of the nearby, cosmological-luminosity GRB 030329 
(M sa 6 x 1O 7 M , similar to that of the SMC) and 
in samples of 0.3 < z < 1 and z ~ 1 GRB hosts 
(jLevesaue et al.ll2010at iCastro Ceron et al.ll2010D . 

The stellar masses of obscured GRB hosts (Figure ITTJf 

27 "Typical" in terms of detected rates — while rarely observed, 
low-luminosity events are more common in a volumetric sense (e.g., 
ISoderberg etal]\200Th . 
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Fig. 10. — Stellar mass versus redshift of GRB hosts and 
field galaxies; symbol conventions are the same as in Figure [7J 
with curves indicating quartile boundaries in the expected distri- 
bution for K < 23 mag galaxies. Dust-obscured GRBs probe a 
much more massive host population than unobscured GRBs, an 
effect which produces the magnitude trend seen in Figure [7J The 
connection between the GRB rate and the SFR appears signifi- 
cantly influenced by mass, especially at z < 1.5 where the relative 
numbers of GRBs in "high-mass" (> 10 10 Mq) versus "low-mass" 
(< 10 10 Mq) galaxies differ by a factor of ~ 5 despite an expected 
ratio of unity. While even deep field surveys are incomplete at 
the low-mass end, very few GRB hosts would have been missed in 
surveys to this depth, and this does not explain our result. Note: 
some additional scatter has been added to the MODS points to 
remove some minor nonphysical discreteness in the mass distribu- 
tion introduced as a result o f the mass-fitting procedure in MODS 
IKaiisawa eTaTl (f20Tol . \20Wi . 

are much larger than those of the unobscured popula- 
tion, a trend manifest in the remarkable fact that all 
dark hosts in our sample are detected in the NIR and 
with IRAC, a reversal of previous conclusions f ocused 
on optically bright bursts ()Le Floc'h et al.l 12006). The 
dark-host population in our sample has a median mass 
of about 5xl0 9 M Q atzwl and 3 x 10 10 M at z « 2, 
roughly an order of magnitude higher than the equiva- 
lent values for unobscured GRBs 28 . We know of only a 
few obscured GRB hosts at z « 3, but comparisons to 
even the maximum mass values derived for the host pop- 
ulatio n of unobscured GRBs at z = 3-4 by lLaskar et al.l 
(|2011h suggests that this difference persists up to at least 
this redshift range. 

While unobscured GRB host galaxies show an obvi- 
ous aversion for the most massive host galaxies identi- 
fied in field surveys (only a single unobscured-GRB host 
has a mass > 2 x 10 10 Mq), dark GRBs show no such 
restrictions and populate both high-mass and low-mass 
galaxies, especially at z > 1.5 where massive hosts are 
quite common. At the same time, it is notable that at 
2«1 we still do not identify any bursts within our sam- 

8 For reference, the stellar mass of the LMC is ~ 3 x 10 9 M 
and that of the Milky Way is ~ 5 X 10 10 M 



Fig. 11. — Average dust attenuation for GRB hosts, as given 
by the results of our SED fitting. Symbol conventions are the 
same as in Figure 1101 with curves indicating quartile boundaries 
of the expected host distribution for K < 23 mag galaxies. Dust- 
obscured GRBs (unlike optically selected GRBs) tend to inhabit 
dust-obscu red galaxies, although the trend is not universal — see 
also Figure [151 

pie above ~ 2 x 10 10 M Q , a mass scale which still con- 
tributes substantially to the SFR density at this epoch. 
Such galaxies are very easy to detect at this redshift (two 
are present in the pre- Swift sample), and there is no rea- 
son we should be missing these objects if they are indeed 
common among GRB hosts. We will return to this ques- 
tion in more detail in Sj7l 

6.3. Average Attenuation 

Given that our sample was selected on the basis of an 
obscured sightline to the host progenitor and that pre- 
vious GRB host samples rarely showed significant host 
extinction, it would be surprising if similar signatures of 
extinction were not manifest in the host galaxy itself. 
In Figure [11] the inferred extinction of each event in our 
sample is plotted against the comparison sample. 

Indeed, the typical dark GRB host is quite dust ob- 
scured. A rest-frame Ay ~ 1 mag is typical, with extinc- 
tions often extending to ~ 2 mag (although not much 
beyond). At the same time, several hosts show no ev- 
idence of extinction at all, so this trend is hardly uni- 
versal. These properties seem to differ from those of the 
unobscured GRBs in our sample, which typically have 
Ay ~ 0.5 mag and rarely exceed Ay > lmag (although 
this could, in part, reflect the redshift differences between 
the populations). 

Most high-z star formation occurs in fairly dust- 
obscured regions. Most of the luminous galaxies that 
dominate cosmic star formation at z = 1-4 are moder- 
ately dust-obscured and visibly dust- re ddened in the UV 
(jMadau et alJll998HMeurer et al.lll999D . with UV atten- 
uations of a few magnitudes. The Ay values observed in 
our dust-obscured hosts appear fairly typical of the sorts 
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Fig. 12.— Dust-corrected UV SFR (as inferred from SED fit- 
ting to the UV/optical/lR observations) for GRB hosts and field 
galaxies. Symbol conventions are the same as in Figure 1101 with 
curves indicating quartile boundaries of the expected host distri- 
bution for K < 23 mag galaxies. The symbol size for the MODS 
survey galaxies is proportional to the total star formation including 
contribution from 24 fim-derived SFRs, although for most galaxies 
this is similar to the UV dust-corrected SFR plotted on the ordi- 
nate. Dust-obscured and optically selected GRBs occupy galaxies 
with a range of SFRs, although rapidly star-forming galaxies are 
much more likely to host obscured GRBs. Several GRBs are hosted 
in probable ULIRGs given their extremely large SFRs. Note that 
there are very few non-dark comparison objects at z > 2 since there 
is no way to determine SFRs from t he one-filter and two -filter pho- 
tometry of lHiorth et~aTI l|20lD and lLaskar et al.l (goTJ). 

of galaxies that dominate the SFR density, whereas the 
hosts of optically selected GRBs appear not to, preferring 
much more optically thin galaxies. 

6.4. Star- Formation Rate 

Inferring the current SFR in a galaxy from photom- 
etry alone is intrinsically more complicated than deter- 
mining most other observable properties because of its 
sensitivity to the star-formation history: a recent shut- 
off (or surge) in the SFR is difficult to recognize, since 
stellar populations continue to produce significant UV 
flux during the first 50-100 Myr. All UV-continuum cal- 
ibrations rely implicitly or explicitly on the assumption 
of constant SFR (or some other prescribed slowly chang- 
ing star- formation history), as we also do here (at least 
during the last 50 Myr; 2J) In addition, the UV bands 
are strongly affected by interstellar dust; for a typical 
luminous z w 2 galaxy, 75% or more of the UV light is 
likely to be absorbed by dust before it is able to escape 
the galaxy. 

In spite of these uncertainties, UV-based estimates of 
SFR are a standard part of the extragalactic toolkit and 
generally have been quite successful at matching SFR 
calibrators at other w a velengths (e . g.jElbaz et al.[l2007t 
Pannella et all 120091: iSalim et al.l 120091: iReddv et all 
2012T ). The submillimeter galaxies (and low-redshift 
ULIRGs) represent an exception to this — in these galax- 
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Fig. 13. — Specific SFR for GRB hosts. Symbol conventions are 
the same as in Figure [Tu1 with curves indicating quartile boundaries 
of the expected host distribution for K < 23 mag galaxies. There 
are no obvious trends visible, and the uncertainty of many of our 
points is large. A significant fraction of both obscured-GRB and 
unobscured-GRB hosts occur in galaxies with very high sSFRs, 
with ages clearly younger than the minimum permitted within the 
MODS modeling. However, this population is not dominant (per- 
haps 20—30% of GRB hosts). A few hosts actually show signifi- 
cantly lower sSFRs, even below l/<H u bble ( re< 3 curve), indicating 
a mature galaxy whose apparent SFR has fallen with time. 

ies most of the star formation occurs in optically thick 
regions; its presence can be inf erred only from the ree mit- 
ted light at long wavelengths ()Chapman et al.ll2005|) . 

Most of the host galaxies in our dust-obscured sam- 
ple show very high SFRs. The median SFR for dark 
GRB hosts at z = 1-2 is ~ 10M Q yr _1 ; for those at 
z > 2 it is ~ GOMQyr^ 1 . This appears significantly 
higher (by a factor of a few) than for the hosts of opti- 
cally selected GRBs, although given the rapid evolution 
of SFR with redshift and the small number of published 
non-dark GRBs at z > 1.5 it is difficult to be as cer- 
tain of this trend as for the other parameters discussed 
previously, since the TOUGH and Spitzer high-z samples 
cannot be used to constrain the SFR even approximately. 
At minimum, however, it is clear that after inclusion of 
dark-GRB hosts, the supposed de arth of very luminous 
galaxies (LIRGs and ULIRGs; see lLe Floc'h et al.ll2005fl 
compared to field samples is at least reduced if not nec- 
essarily alleviated completedly. 

6.5. Age and Specific Star- Formation Rate 

We are not actually able to constrain the age directly 
for most of our sample beyond distinguishing crudely be- 
tween very young galaxies with a weak or absent Balmer 
break (and therefore a very young age) and galaxies with 
a strong Balmer break (and therefore an age of a few hun- 
dred Myr or more) in the case of hosts that are not too 
obscured — except for a few objects, we lack the high-S/N 
photometry around the Balmer break to precisely mea- 
sure it. However, since we can measure the mass and 
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SFR, this constrains the age indirectly in terms of the 
specific SFR (sSFR = SFR/M) and its inverse, which 
is effectively a characteristic timescale for forming the 
galaxy if star formation proceeded at its current rate 
since formation. 

The sSFRs of our sample are plotted in Figure [131 Val- 
ues of sSFR r; 1 Gyr -1 (equivalently, formation times of 
about 1 Gyr) are typical for both dark and non-dark 
hosts, seemingly regardless of redshift. However, there 
are significant excursions in both directions: a few dark 
GRBs have sSFRs about one order of magnitude above 
this, indicating very young (dominant) stellar popula- 
tions. Given the short timescales involved (10 7 -10 8 yr), 
these galaxies clearly are in the throes of a short-lived 
event such as a galaxy merger. On the other side of the 
plot, older galaxies seem to produce dark GRBs almost 
exclusively: with perhaps one or two marginal excep- 
tions, the only hosts with inverse sSFRs equivalent to 
the Hubble time (indicating long-term quiescent evolu- 
tion over many Gyr or a declining sSFR) are those of 
dark GRBs. 

Other than this, no obvious trends are visible. Several 
hosts reach above the limit of maximum sSFR of the 
MOI RCS Deep Survey catalog of iKaiisawa et all (|2010l 
|2011|) . although this is not necessarily a physical state- 
ment as the assumed star-formation history in the MODS 
modeling did not allow for younger stellar populations 
than 50 Myr. Much more insight can be gleaned from 
examination of the sSFR in relation to other parameters, 
which we will return to in the discussion. 

7. DISCUSSION 

Previous, smaller-scale studies of the hosts of dark 
GRBs have led to some disagreement regarding the na- 
ture of this group of galaxies relative to the broader 
GRB host population. In many cases, the detection 
of a relatively blue, unobscured host has led to sugges- 
tions that even the inclusion of dark GRBs is unlikely 
to change the overall GRB host population s ubstan- 
tially (e.g., iGorosabel et al.1 l2003t [Tanvir et all l2008bb 
iPerlev et alJ l2009). But other studies have found very 
different results, with the detection of luminous, red host 
galaxies possibly necessitating a wholesale revision of 
our notion of w hat types of galax ies h ost GRBs (e.g., 
Hashimoto et alJl201oYlHunt et alj|2011t ISvensson et al.1 
2012tlKriihler et al.ll2012bURossi et al.ll2012h . 

Given the above results, the basic reason for this con- 
tradiction is now apparent: the population itself is intrin- 
sically very diverse, and small numbers of objects neces- 
sarily show only a subset of this diversity. In fact, the 
host galaxies of dark GRBs cover many orders of magni- 
tude in almost every property, ranging across most of the 
parameter space encountered among field-selected high-z 
star-forming galaxies detected in current surveys. Com- 
pounding this, the properties of both obscured and un- 
obscured GRB hosts appear to evolve substantially with 
redshift, confusing any comparisons unless redshift dif- 
ferences in the populations are taken into account. 

7.1. Do GRBs Follow Star Formation? 

Given this diversity — and the fact that the dark GRB 
host population seems to represent an almost perfect 
complement to the host population of low- Ay, optically 
bright GRBs (Figure IT4|) — it is tempting to rush to the 



conclusion that GRBs follow star formation in a uniform 
way after all. Broadly speaking, every type of high- 
z star-forming galaxy appears capable of producing a 
GRB: the full range of masses, extinctions, SFRs, and 
other parameters now seems to be represented among 
the GRB host population. 

However, the question is not just whether all types 
of galaxies can produce GRBs but whether they do so 
at the rate that we expect. Only if GRB formation is 
strictly limited to a small region of parameter space by a 
critical parameter (e.g., by a hard metallicity cutoff at a 
specific value) will GRBs avoid certain regions of galaxy 
parameter space entirely. The question is one of number 
statistics, matching the GRB rate as a function of various 
parameters (redshift, mass, luminosity, etc.) to the SFR. 

To illustrate the point more clearly, in Figure [14] we 
plot the mass of objects in our sample against the UV- 
inferred optical SFR (which is currently all we are able 
to measure without longer- wavelength data) . The plot is 
divided into four redshift bins: z < 0.5, 0.5 < z < 1.5, 
1.5 < z < 2.5, and 2.5 < z < 3.5 respectively (SDSS 
is used as the comparison sample at z < 0.5 instead of 
M ODS, using the DR7 M PA-JHU value-added catalog 
of lBrinchmann et al.ll2004l ) The distinction between the 
dark and non-dark hosts is easily apparent, and the in- 
clusion of dark GRBs greatly improves the consistency of 
the host population with the overall star-forming galaxy 
population probed by the MOIRCS Deep Survey. 

Since GRBs are star-formation selected, they should 
not occupy the same space filled by the field galaxies in 
the diagram uniformly but should prefer galaxies with 
higher SFR in direct proportion to their SFR. We have 
represented this fact visually on our plots by scaling the 
area of the MOIRCS galaxies in proportion to their total 
SFR — loosely speaking, on our diagrams we expect the 
GRB distribution to follow the density of gray "ink." In 
fact, even with dust-obscured bursts included, this is not 
how the GRB host population actually behaves: at least 
at lower redshifts (in particular the 0.5 < z < 1.5 bin), 
the top end of the diagram seems to remain underpopu- 
lated, with most events crowding toward the bottom and 
the left, indicating low masses and SFRs. While the dust- 
obscured bursts are more luminous and massive, they 
are not sufficiently different to reverse this trend. Nor 
are they sufficiently common: events of the type ana- 
lyzed here contribute only about 15% of the population 
(Figure [T]), much less than the ratio shown in the figures 
above z > 1. 

7.1.1. As a Function of NIR Magnitude 

The above argument is not quantitative and, impor- 
tantly, it ignores the incompleteness of the comparison 
survey. By definition, flux-limited galaxy field surveys 
only probe star formation above a certain limiting flux 
level, while GRB-targeted surveys probe all star forma- 
tion including that in galaxies below the flux level (which 
will show up as nondetections if observed to the same 
depth as the field survey). This will produce a corre- 
sponding difference in the overall galaxy populations that 
will complicate detailed comparisons of the two samples 
if not corrected. 

This problem can be avoided easily by simply mimick- 
ing the cut set by the survey depth threshold on the GRB 
host sample. For the MOIRCS Deep Survey, the selection 
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Fig. 14. — SFR versus mass for obscured and unobscured GRB host galaxies compared to galaxies from SDSS and the MODS field 
survey at a variety of redshift ranges: z < 0.5, j fs 1, z « 2, and z fa 3. (An SDSS slice of 0.1 < z < 0.15 is used as the comparison sample 
at z < 0.5; data are taken from the MPA-JHU catalog.) We also show mass- SFR correlations roughly corresponding to each redshift range 
USalim et al. 2007; El baz et all 12007: Daddi ct al. 2007; Magdis ct al. 2010s). Despite the large sample of dark GRB hosts added to the 
population, we continue to see a notable deficiency of high-Af, high-SFR galaxies at Z 1 (if not necessarily at higher redshifts). Symbol 
conventions are the same as in previous plots. 



criterion was a magnitude cut at K < 23.1 Vega mag (the 
5cr threshold for their "wide" field). This limit is signif- 
icantly fainter than the flux of any of our dust-obscured 
hosts, suggesting that all of those galaxies would be re- 
covered in a MODS-likc field survey. Many non-dark 
GRBs, however, are in galaxies that are undetected in K, 
and these upper limits do not actually reach K w 23 mag 
in any case (which is a very deep limit, difficult to achieve 
in practice). However, we can calculate their expected 
-ftT-band fluxes approximately via synthetic photometry 
of our best-fit model (in essence, extrapolating the flux 
measurements in the more-sensitive bluer filters). Based 
on this procedure, we find that all except for two of the 
GRB host galaxies in the (mostly lower-redshift) pre- 
Swift comparison sample would also have been recovered 
in MODS, indicating that — as advertised — this very deep 
survey reaches flux levels that include the large major- 
ity of the star formation probed by GRBs at z < 1.5. 



Consequently, field-survey incompleteness has almost no 
effect on the comparison between our results and this 
survey 29 . Its effects can be eliminated completely by re- 
moving these two galaxies from consideration. 

We illustrate this point explicitly in Figure [71 The sur- 
vey threshold is shown as a dashed horizontal line; nearly 
all GRB hosts we consider lie above this level. This figure 
also includes a curve showing the median SFR-weighted 
galaxy magnitude of the comparison sample as a function 
of redshift (i.e., for a given redshift, we have calculated 
the magnitude at which 50% of star formation probed by 

29 This is, in part, due to the fact that the pre- Swift compar- 
ison sample contains very few events at z > 2. A higher-redshift 
comparison sample would almost certainly contain a much larger 
number of g alaxies below the M ODS threshold — as evidenced by 
the work of Laskar ct al. (2011), which we will consider in more 
detail in the next section. 
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the survey 30 occurs in brighter galaxies and 50% occurs 
in fainter galaxies). Dotted lines show the 1st and 3rd 
quartile boundaries (i.e., the 25th and 75th percentiles). 
(To smooth out statistical variations and effects of cos- 
mic variance and in this small-angle survey, we have fit a 
polynomial over the full redshift range for each line.) If 
GRB hosts are to represent uniform star-formation trac- 
ers, then the numbers of events should divide approxi- 
mately evenly between the four lines, except below the 
bottom horizontal line at K > 23 — numbers in galaxies 
fainter than this limit cannot be predicted since the field 
survey does not probe beyond this depth. (To empha- 
size this, we plot host galaxies fainter than this level with 
hollow points.) 

The comparison is most easily made at z ~ 1, where 
there is a large pre- Swift GRB comparison sample and 
where the incompleteness limit of MODS (and typical 
observed upper limits on faint GRB host galaxies) are 
far below the magnitudes of the galaxies at which most 
star formation is actually occurring. We focus specifi- 
cally on the redshift range of 0.6 < z < 1.4, which con- 
tains 31 events in total between our sample, the pve-Swift 
sample, and TOUGH; 7 of these are "dark," a fraction 
that is probably comparable to the actual co ntribution of 
obscured bursts to the o verall population (|Cenko et al.l 
120091 iPerlev et all 120091 : iGreiner et aTll2011ft . Of these 
events, only three are brighter than the expected me- 
dian line and none are brighter than the 75th-percentile 
line! All other GRBs (except for the five nondetected 
TOUGH hosts) land in the bottom quartile but above the 
magnitude limit of the MODS survey, meaning that their 
overabundance is not simply due to GRBs probing star 
formation below the survey limit — GRBs genuinely seem 
to be overabundant in faint galaxies and underabundant 
in bright ones at this redshift. The statistical likelihood 
of producing an empty top quartile in a sample of this 
size (ignoring the TOUGH nondetections) is 0.05%. A 
K-S test between the observed if-magnitude distribution 
for detected galaxies at 0.6 < z < 1.4 relative to the ex- 
pected distribution for a SFR-weighted sample with the 
same redshifts (for K < 23 galaxies) gives a probabil- 
ity of 0.003% that the observed distribution is consistent 
with random chance. 

This trend weakens significantly toward higher red- 
shifts: at z > 1.5 significant numbers of GRB hosts 
do populate the high-luminosity region of the plot and 
agreement with the assumption of an unbiased SFR- 
tracing population improves. It is even possible that the 
"bias" disappears entirely — host galaxies detected in K- 
band in the uniform TOUGH survey between 1.5 < z < 
2.5 (the small circles in Figure excluding upper lim- 
its) are distributed almost uniformly between the four 
SFR-weighted quartiles, seemingly quite consistent with 
an unbiased SFR-tracing population. On the other hand, 
most of the nondetections in the TOUGH survey are 1-2 
magnitudes short of the MODS limit and a large popula- 
tion of faint (but MODS-detectable) hosts may have been 
missed by that effort. In Figure [8] we compare instead 
against a much deeper survey, the HST Snapshot project 

30 We emphasize that these median and quartile boundaries only 
probe star formation in galaxies above the K > 23.1 mag survey 
threshold. Since we have discarded any GRBs occurring within 
fainter galaxies from consideration in this analysis, star formation 
in galaxies fainter than this limit does not concern us. 



(Tibbetts-Harlow et al. in prep.), a randomly-selected 
sample of the hosts of optically-bright GRBs observed 
to very deep limits with WFC3-IR on HST. While the 
number of events in the z ~ 2 redshift range is still rela- 
tively small, five out of seven hosts at F1QQW > 23.5 AB 
mag (roughly corresponding to the typical upper limit 
of K > 21.5 Vega mag in TOUGH) are in fact far be- 
low that threshold and would not have been recovered 
in MODS, suggesting that most of the TOUGH nonde- 
tections might also fall below the MODS survey limit 
and therefore not affect the comparison between the de- 
tected hosts and the MODS survey. However, both com- 
parison surveys are small, and it is quite possible that a 
larger sample would reveal a clearer trend towards fainter 
galaxies at z ~ 2 similar to (if weaker than) what is seen 
at z ~ 1. (It is also worth noting that in terms of color, 
even the most luminous hosts from the uniform TOUGH 
survey appear unusually blue relative to expectations for 
az~2 SFR-selected population; Figure [9l) 

We emphasize that this result is not dependent in any 
significant way on the details of our SED modeling or fit- 
ting procedures — it is a simple matter of number counts 
as a function of K magnitude. Our modeling is used only 
to extrapolate the SED for the faintest galaxies to con- 
firm that they are not so faint that they would not be 
recovered in deep field surveys. Similarly, it is not partic- 
ularly sensitive to the use of the pre- Swift host compari- 
son sample — it is noteworthy that even the population of 
dust-obscured GRB hosts alone, despite being relatively 
massive compared to unobscured GRBs, does on its own 
still trend toward subluminous galaxies at z w 1. 

Our conclusion is of course dependent to some degree 
on the modeling and cosmic volume of the MODS survey 
itself. However, the strength of the discrepancy that we 
observe would require very large errors to explain. In 
particular, the actual rate of star formation in faint (sub- 
25th percentile) z w 1 galaxies would need to be a factor 
of 2-4 higher than claimed relative to the SFR density in 
bright galaxies. Dust extinction in faint galaxies is not a 
viable explanation for why this hypothetical extra SFR 
was missed from MODS, since if this were the case the 
number of dust-obscured GRBs would be much larger 
than the ~ 20% value that is currently reported. Cosmic 
variance, in turn, is unlikely to significantly affect even 
a very narrow-angle slice covering such a wide redshift 
range (~ 2 comoving Gpc between z 0.5 and z « 1.5). 

Finally, we consider the possibility that our GRB sam- 
ple is somehow still nonrepresentative of GRB host galax- 
ies in a significant way. Clearly our "comparison" sample 
is quite inadequate at z w 2 and beyond, but questions 
could also be raised about the biases involved in the re- 
porting of the pre-Swift z 1 comparison sample. Nev- 
ertheless, it is extremely difficult to imagine a situation 
in which our z « 1 sample would be biased to such an 
extent to reproduce the trends observed. For example, 
we have illustrated as an example that out of 31 known 
GRBs we have considered in the range 0.6 < z < 1.4, we 
identify no events in hosts brighter than the 75th SFR 
percentile curve. Hosts of that brightness (K s=s 19 mag 
at z = 1) are easy to find, even with relatively small 
(1-2 m) telescopes. Yet not a single host with K < 19 
mag and z > 0.9 has been reported in the entire pre- 
Swift era or within the unbiased TOUGH survey (which 
includes 10 events in the 0.6 < z < 1.4 redshift range, 
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none of which are brighter even than the expected 50th 
percentile in K magnitude) — indeed, to our knowledge, 
no such hosts have been reported in the entire Swift era. 
Even if there are somewhat more "dark" GRB hosts than 
we realize at z « 1, even these events are just as likely 
to occupy faint galaxies at this redshift than bright ones, 
and so would only slightly improve the consistency be- 
tween expectations and observation. 

The conclusion that GRBs do not trace cosmic star 
formation (at least at z < 1.5) therefore appears very 
difficult to avoid. 

7.1.2. As a Function of Stellar Mass 

While the analysis above in terms of observed mag- 
nitude provides strong evidence of a difference between 
the GRB rate and the SFR, it offers no direct insight 
into its physical origin. However, the same argument 
presented above in terms of the observed K magnitude 
can be recast in terms of physical properties quite easily. 
While addressing the relative completeness of the pop- 
ulations is less straightforward, the fact that over 95% 
of our host population would in fact be recovered to the 
MODS limit implies that this technique is still powerful 
when employed on a range of variables. To avoid the 
risk of incompleteness, we simply discard the two events 
which would have fallen below the MODS threshold (they 
are plotted as hollow points). 

The solid and dashed curves in Figure [TU] correspond 
to the same concepts used in the magnitude-based analy- 
sis, indicating 25th, 50th, and 75th percentile of the SFR 
density as a function of mass at each redshift. Unsurpris- 
ingly (since the NIR luminosity is a close tracer of stellar 
mass) this diagram shows a very similar result as the 
if-band analysis, with most GRB hosts (including many 
dust-obscured GRBs) falling below the 25th percentile 
line at z ~ 1. 

It is worth noting that while the conclusion of 
iKajisawa et al.l 1)2011 ') that the cosmic SFR is domi- 



nated by moderate- to high-mass galaxies (> 10 10 M Q ) 
at z ~ 1 has been supported by o ther recent work 
([Zheng et al.ll2007t ISantini et a,l.ll2009h , some other stud- 
ies (in particularly optically selected samples) do suggest 
a large r role for star forma tion at lower masses: for ex - 
ample, Uuneau et all (|2005| ) and iMobasher et al.l (|2009h 
suggest that the median mass is somewhere near or be- 
low 10 10 M Q at z w 1-2. However, the much shallower 
de pths of these survey s, K < 20.6 mag (Vega) in the case 
o f lluneau et all 12003 and i < 25 mag (SDSS/AB) for 



MobashereLaL|2009, makes them become incomplete at 
around 10 10 M Q , making it difficult to actually evaluate 
the amount of star formation occurring below these levels 
where most of our sample lurks. Purely optically derived 
SFRs also miss contributions from very dusty galaxies 
which, while not necessarily dominant, are at minimum 
an important contributor to the SFR density at higher 
redshifts. 

Mass, therefore, appears to be a primary driver of the 
deviation between the expected and observed character- 
istics of the GRB host population. Assuming the balance 
of "unobscured" and "dark" GRBs in our combined sam- 
ple at z w 1 is close to the intrinsic one, the difference 
in GRB rate per unit star formation varies by a factor of 
at least 5 and possibly by 10 or more between 10 11 M 
galaxies and 10 9 M Q galaxies. This observation is eas- 



ily interpreted under the current (if still-controversial) 
paradigm that the GRB progenitor is strongly metal- 
averse, since the mass and average metallicity of a galaxy 
are strongly correlated. Nevertheless, mass correlates 
with a number of other observables as well, and so it is 
useful to check to determine what other bulk properties 
may contribute to the observations. 

7.1.3. As a Function of Luminosity and SFR 

It is natural to ask how the cosmic SFR density 
is divided among the individual SFRs of the galaxies 
that contribute to it, including the relative balance of 
the large numbers of individually insignificant low-SFR 
galaxies versus spectacular but rare systems such as sub- 
millimeter galaxies, ULIRGs, and the most luminous 
LBGs. The actual answer to this question has remained 
controversial over the years and (often) highly dependent 
on the wavelength at which the study was conducted, al- 
though the most recent studies have been converging on 
a contribution of ~ 20% from highly luminous and ob- 
scured super-star- forming galaxies (ULIRGs and SMGs) , 
with the remaining fraction about evenly distributed be- 
tween moderately luminous galaxies (LIR Gs and lumi- 
nous L BGs) and lower-lumino si ty systems (iReddv et al.1 
20081 iReddv fc Steidell l2009t iRodighiero et al.l ]201oF 



Murphy et al.l 120111: iMagnelli et al.l 1201 ll ) . The impor- 
tance of luminous galaxies peaks at z = 2-3 and falls 
rapidly toward lowe r and (probably) higher redshifts 
(jBouwens et al.ll2009D . 

We are currently restricted to the dust-corrected UV as 
a probe of the SFR, which is insensitive to star formation 
in high-optical-depth regions. Long-wavelength observa- 
tions are being pursued and will be the basis of future 
work, but even without those constraints the problem 
can be framed in terms amenable to the data at hand 
using the same basic technique as in previous sections. 
In Figure[l2]we calculate the median and 25th/75th per- 
centile boundaries for all star formation in K < 23 mag 
galaxies as a function of the UV-inferred, dust-corrected 
SFR, which corresponds to the same measure probed by 
our sample. Despite reports of a GRB sample deficie nt in 
luminous galaxies and LIRGs (|Le Floc'h et al.ll2006l ). we 
actually observe reasonably good agreement in this prop- 
erty with a significant number of events occupying the 
most luminous quartile at all redshifts, and even with- 
out direct long-wavelength constraints we predict (from 
the UV luminosities and mean dust attenuations) that 
most of our dust-obscured events are LIRGs and a few 
are (borderline) ULIRGs. At z rs 1 where we previously 
observed a clear trend among GRB hosts toward low- 
mass galaxies, we infer only a slight preference toward 
lower-luminosity galaxies, and this deviation is consistent 
with statistical variations. (Again, at z > 1.5 the non- 
dark comparison sample is inadequate at present to claim 
that our overall sample is representative of all GRBs.) 

Given that samples of high-z gala xies show a corre- 
lation between mass and SFR (e.g., iDaddi et alJl2007t 
lElbaz et al.l l2007h . the fact that the trend observed in 
stellar mass does not carry over to SFR is somewhat sur- 
prising, since even if SFR were physically unconnected to 
the GRB rate we would somewhat expect a tendency of 
low-mass galaxies to couple to low-SFR galaxies for this 
reason. It is possible that a larger sample would recover 
this trend. Alternatively, it could be actively counter- 
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acted by a third parameter — in particular, a tendency 
for GRBs to prefer high-sSFR galaxies, a possibility we 
explore in the next section. 

7.1.4. As a Function of Specific Star- Formation Rate 



As we noted in A6.51 a remarkable number of our host 
galaxies show very high sSFRs, evidenced directly by 
young SEDs with little or no hint of a Balmer break, 
indicating average stellar ages younger than 100 Myr or 
(equivalently) sSFRs of > lOGyr" 1 . About 20% of our 
sample lies above this threshold. This tendency for GRBs 
to occur in young, high-sSFR galaxie s has been noticed 
previously among non-dark GRBs ((Ch ristensen et al.l 
12004 iMannucci et alJl201lHKocevski k WestJl2011l ): the 
inclusion of our sample of dark GRBs extends the trend 
up to higher redshifts (and stellar masses), as shown in 
Figure [HI 

It is important to note that the association of GRBs 
with very young galaxies is far from exclusive — the re- 
maining 80% of the population we examined occurs in 
galaxies with a quite "normal" population ages for a 
SFR-selected population. Nevertheless, it is interesting 
that a very significant fraction of GRBs occurs in ex- 
tremely young galaxies whose contributio n to overall cos- 
mic st ar formation is almost negligible: Kaiisaw a et al.1 
(|2011[ ) find only a few percent of all SFR density oc- 
curring in > lOGyr -1 systems at z < 2, most of 
which comes from galaxies with 10 10 < M < 10 10 ' 5 Mq. 
The fraction rise s with redshift, but is still low even 
at z ^ 2 (10%; iRodighiero et al.l l20lU although see 
lYoshikawa et al.ll201Ct) 

This is not just a reflection of the mass difference dis- 
cuss ed earlier (^7.1.21)— the mass-sSFR relation is nearly 
flat jElbaz et alJl2007t iDaddi et al J2007L IWhitaker et all 
12012ft . and so (at a given redshift) the apparent skew 
toward lower mass should not produce any comparable 
differences in sSFR. Rather, for a galaxy of a given mass 
at a given redshift, the GRB rate in high-sSFR galaxies 
seems to be substantially elevated above what is expected 
given the SFR. The origin of this trend is not clear; while 
such an effect is predicted by th e three-paramete r mass- 
sSFR -metallicity correlation of iMannucci et al.l (|2010L 
1201 lh , in practice the sSFR-metallicity portion of the cor- 
relation is fairly weak and unlike ly to produce the degree 
of the correlation observed here (jKocevski k Westll2011l ) 
unless the actual metallicity dependence of the GRB rate 
is quite strong (and possibly continuous across the en- 
tire mass range — note that the apparent preference for 
high-sSFR systems is seen across the entire mass range 
of star- forming galaxies, as can be seen in Figure III)) . 
More work will be necessary to determine if this sort of 
metallicity effect would be sufficient, or if dependence 
on another parameter more directly correlated to spe- 
cific star-formation (such as the UV radiation field of 
the galaxy) would need to be introduced to explain the 
apparently elevated rate in these systems. 

7.2. Correlating Afterglow and Host Extinction — The 
Distribution of Host Dust 

While GRBs may not be unbiased statistical tracers of 
star formation, they are still able to address questions 
regarding the detailed internal characteristics of high-z 
host galaxies that cannot be answered by more tradi- 
tional field-survey techniques at all. 
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Fig. 15. — Compar ison of extinction measured along the af- 
terglow sightline (from Kann ct al. 2006 or this work) against the 
average extinction measured in the host from our SED-fitting pro- 
cedure (which generally indicates the average attenuation toward 
young stars in the optically thin portions of the galaxy). The dot- 
ted line indicates equality of the two values. Most of the points arc 
fairly consistent (within a factor of 2-3) with this line, indicating 
that the ISM in these systems is relatively homogeneous and that 
dustier galaxies do indeed tend to produce dustier bursts. A mi- 
nority of systems do show significant excursions, in all directions: 
effectively unobscured GRBs in moderately dusty galaxies exist, as 
do highly obscured GRBs in blue, unobscured galaxies. Such exam- 
ples correspond to more heterogeneous galaxies where the sightline 
extinction to an individual star can deviate greatly from the av- 
erage. Nearly all of the most conspicuously heterogeneous objects 
are of intermediate mass. 

In Figure [TS] we plot the extinction as measured from 
the afterglow (determined either by our own preliminary 
analysis in j^l fo r dark-burst hosts or by the study of 
iKann et "al"1l2006l for pre- Swift hosts) against the inferred 
average Ay of the galaxy's starlight (in optically thin re- 
gions) from our host SED fitting. If every star in the 
galaxy were obscured by the same amount and type of 
extinction (a homogeneous screen), the points would all 
fall along the diagonal line plotted. Excursions from this 
line (in either direction) indicate heterogeneity in the 
dust distribution with respect to the stars, while a pref- 
erential excursion would indicate a tendency for GRBs 
to inhabit either the dustier or less dusty parts of the 
host (relative to the typical optically-thin region). 

Most points in the diagram are in fact fairly close to 
the diagonal line (within a factor of ~ 2), indicating 
that in most galaxies the extinction is constistent with 
a fairly homogeneous distribution that affects most stars 
in a similar way. This result confirms that the simpli- 
fied picture used to treat extinction in populations of 
high-redshift galaxies (namely, the treatment of extinc- 
tion and reddening of galaxies as a simple screen with 
a roughly universal attenuation curve) is at least rea- 
sonable in many objects, and suggests that, when the 
amount of obscuration seen along a GRB sightline is 
low to moderate (0.2-2 mag), this dust is probably dom- 
inated by the diffuse ISM of the galaxy (since most other 
stars are attenuated by a similar magnitude ). Given pre- 
vious claims (by our group and others; e.g.. iPerlev et all 
2009) that extinction may merely be a local effect, this 
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correlation is remarkable. 

On the other hand, this result is far from universal — a 
large number of events deviate far from the basic expec- 
tation of a linear relation in Figure [T5l In particular, the 
GRB sightline extinction can be much, much larger than 
the extinction obscuring a "typical" star in the optically 
thin portion of the galaxy — indicating that the GRB lies 
inside or behind a region of optically thick dust that is 
not representative of the diffuse ISM filling the rest of the 
galaxy. Typically, galaxies hosting these extremely ob- 
scured GRBs are already moderately dusty overall (the 
GRB sightline is just dustier than average), but in two 
cases there is clear evidence of a very heavily obscured 
GRB occurring in a galaxy with basically no evidence of 
obscuration at all. Such events are rare and appear to be 
outliers from the rest of the population, but there is little 
question of their existence. (One event, GRB 061222A, 
is practically a prototype of the heavily obscured Swift 
GRB: NIR-only detections, very red H — K color, large 
Nn, unambiguous galaxy association, host spectroscopic 
redshift.) 

There are several possible interpretations of this obser- 
vation, as follows. 

(1) Extremely heterogeneous ISM. One possibility is 
that these (relatively uncommon) associations of ob- 
scured bursts in blue galaxies indicate random align- 
ments between the GRB progenitor and an unassociated, 
dense molecular cloud in an otherwise mostly transpar- 
ent galaxy. Something like this probably happened with 
GRB 080607 — spectra of the optical afterglow of this 
heavily obscured GRB indicate that the burst sightline 
penetrated a discrete molecular cloud along its sightline 
at a distance of a few hundred pc from the burst it- 
self (as evidenced from the excitation of fine-structure 
lines a nd nondestruction of molecular lines in th e spec- 
trum; iProchaska et al.1 l2009t iSheffer et all l2009t ) . The 
discrepancy between host and afterglow Ay for this GRB 
(1.1 mag vs. 3.3 mag, respectively) is not as large as for 
some other examples, but could indicate that this sce- 
nario is widespread. 

(2) Local dust. The most recent star formation in 
most galaxies is obsc ured significantly mor e heavily than 
slightly older stars ([Chariot fc Falll |200Q[ ). since newly 
formed stars often remain embedded in the (dusty, opti- 
cally thick) molecular cloud that formed them for a few 
Myr after formation, whereas the older stars that domi- 
nate the continuum UV flux have time to disperse these 
clouds entirely. A sufficiently massive, short-lived pro- 
genitor star may explode on a short enough timescale 
that the surrounding gas and dust are not yet cleared. 
Since only a few GRBs are heavily obscured this effect 
cannot be very significant in most cases (either the op- 
tically thick covering fraction is quite small or the GRB 
event is able to disperse the cloud on its own via its 
X-ray and /or UV emission: iWaxman fc Draind 120001: 
iFruchter et "all 120011: iPerna et all I2003D . but under the 
right circumstances such local dust could produce a large 
obscuring column. 

(3) Nuclear super- starbursts. The extinction inferred 
from the UV-continuum only probes optically thin re- 
gions of the galaxy. In fact, if a small, extremely heavily 
obscured portion of the galaxy were the site of most of its 
star formation (as in local ULIRGs such as Arp 220 and 
in submillimeter galaxies), then the Ay and SFR deter- 



mined from the optically thin regions would be entirely 
inappropriate for the galaxy overall. Deep submillime- 
ter or radio observations (which will be presented in a 
companion paper) would be necessary to check this hy- 
pothesis. 

Excursion in the other direction is also present — while 
the typical optically bright GRB has a blue, low-Ay host, 
the non-dark host population does contain a few moder- 
ately obscured galaxies as well. This could simply be a 
reflection of geometry (since stars are embedded inside 
the diffuse interstellar dust and not actually behind a sin- 
gle screen, some should always be much less obscured and 
others somewhat more obscured), though if GRBs pre- 
ferred the outer, less-obscured regions of their hosts (e.g., 
due to an aversion of the more metal- rich center), such 
a trend would be greatly amplified. The submillimeter 
hosts of optically bright GRBs menti oned previously — in 
particular GRBs 0418 and 010222 (iBerger et aLl l2003l: 
iTanvir et all l2004t iMichalowski et all l2008[ )— probably 
represent much more extreme and unambiguous exam- 
ples of this phenomenon. 

The situation is complex, and all of these effects may 
operate, affecting some bursts and not others. More 
observations — in particular high-resolution observations 
of both dark GRBs and their host galaxies, and after- 
glow spectroscopy of fortuitous events like GRB 080607 
luminous enough to be observable optically through their 
dust screens — will be necessary to say more about this 
issue. Nevertheless, although a high-extinction afterglow 
sightline does not guarantee a dusty host and vice versa, 
the fact that the two properties are correlated does imply 
that the origin of most of the dust that is observed for 
at least moderately obscured GRBs (if not necessarily 
extremely obscured GRBs) is ordinary dust in the ISM; 
the darkness of a GRB is not purely geometric. 

7.3. Implications of the Mass-Obscuration Correlation 

Stellar mass is well known to correlate with the average 
attenuation of a star-forming galaxy within field-selected 
samp l es (e.g.. [Brinchmann et all [20041 : iPannella et all 
l2009t iGarn fc Bestl l2010f ) . although there is significant 
scatter in the relation and these variables also correlate 
with many other properties, such as the SFR and metal- 
licity. 

Similar trends are, unsurprisingly, observed in our sam- 
ple. More notable, however, is that the correlation be- 
tween mass and obscuration remains extremely strong 
when comparing host mass to afterglow obscuration, in 
spite of the large number of outliers in the host vs. af- 
terglow comparison of Figure [T5l To emphasize this, we 
have color-coded the points based on the mass of the 
host galaxy, with the most massive hosts (> 4 x 10 10 M ) 
in red, the least massive (< 2 x 10 9 M Q ) in green, and 
intermediate-mass objects in yellow. The mass of the 
host galaxy is clearly a strong predictor of both measures 
of obscuration, but the strong separation between the 
sightlinc-cxtinction observed in the lowest-mass galaxies 
and the highest-mass galaxies is particularly notewor- 
thy. While some caution should be used in this com- 
parison given that the effects of redshift evolution are 
not considered, we find that the least-massive galax- 
ies host almost universally unobscured GRBs while the 
most massive galaxies host almost universally obscured 
GRBs. Intermediate-mass galaxies, in contrast, show 
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much greater diversity and represent all of the strong 
outliers on this plot. No other property we have exam- 
ined (SFR, sSFR, and the Ay of the galaxy itself) shows 
such an obvious correlation at the extreme ends. 

Why would mass be the primary driver of afterglow 
obscuration, more important than even the mean extinc- 
tion in the galaxy itself? Presumably, mass is correlated 
not just with the amount of dust but also with the degree 
of homogeneity or heterogeneity of that dust. Low-mass 
galaxies must be low in dust throughout — there are very 
few sightlines through such objects that produce appre- 
ciable extinction, and so only the occasional rare align- 
ment piercing an unassociated molecular cloud produces 
substantial absorption. On the other end of the scale, 
high- mass (M > 2 x 10 10 M Q ) star- forming galaxies offer 
very few escape routes for a GRB beam to exit the host 
without encountering at least some significant amount 
of dust, suggesting that a significant amount of obscur- 
ing dust is distributed relatively homogeneously within 
the galaxy (as also evidenced by its heavily reddened 
starlight). Of course, such galaxies may also have their 
own localized clouds (or large nuclear starbursts) that 
may increase the attenuation of certain sightlines still 
further. 

This class of intermediate galaxies represents the bulk 
of cosmic SFR (though not, as we have discussed, the 
bulk of GRBs except at z > 1.5) and includes a wide 
diversity of both host colors and afterglow colors, in- 
dicating a group of galaxies that is heterogeneous as a 
population and (often) also internally heterogeneous in 
terms of the distribution of dust and stars within its ISM. 
On the other hand, all of the most massive star-forming 
galaxies in this sample obscure their entire stellar pop- 
ulation (including the individual star that produced the 
GRB), while the least massive galaxies rarely obscure 
any of their stars. 

8. CONCLUSIONS 

We have compiled by far the largest sample of dust- 
obscured GRB host results to date, targeting a total of 23 
GRBs with afterglow-inferred host extinction columns of 
Ay > 1 mag (corresponding to the dustiest ~ 25% of all 
GRBs) with intensive follow-up observations at optical- 
NIR wavelengths both on the ground and in space. This 
sample is compared against a compilation of optically 
selected, well-studied pie-Swift GRB hosts as well as a 
recent deep, narrow-field galaxy survey. 

We find large, unambiguous differences between the 
populations. Dust-obscured GRB hosts are more massive 
than the hosts of optically bright GRBs within the same 
redshift range by about an order of magnitude, and show 
similar tendencies toward higher luminosities and SFRs, 
redder colors and higher dust extinctions, and older stel- 
lar populations on average. 

Despite the revision upward in the average mass and 
luminosity of the GRB host population that these results 
imply, dark GRB hosts are neither sufficiently frequent 
nor massive enough to bring the overall population into 
line with the expectations of a purely SFR-tracing pop- 
ulation at z w 1. The average mass of a GRB host is 
about one order of magnitude lower at this redshift than 
would be predicted from the population of galaxies in- 
ferred from deep multiwavelength field surveys. Consis- 
tency improves significantly at z > 1.5; at these earlier 



epochs massive GRB hosts become fairly common, ap- 
proaching the frequency that would be expected for an 
SFR-tracing population at these redshifts. 

Among the parameters we analyzed, the primary fac- 
tor affecting a galaxy's GRB rate per unit star formation 
is mass — GRBs avoid high-mass galaxies and prefer very 
low-mass ones. Without additional data we can only 
speculate about the cause of these biases at present, al- 
though because mass and metallicity a re well correlated 
(jTremonti et all l2004t lErb et all [20061 ) it seems consis- 
tent with the general idea that metallicity is a funda- 
mental driver of the GRB rate. If so, metallicity is likely 
to act as a gradual modifier of the GRB rate, rather 
than impose a strict cutoff at a specific v alue — for exam- 
ple, w e can strictly rule out the model of iKocevski et all 
(2009), who predict that this metallicity effect should 
impose an equivalent mass cutoff of 1O 9 M at z = 1. 
In reality, of course, a metallicity cut would be some- 
what softened by dispersion in the mass-metallicity re- 
lation as we ll as chemic a l inhomogeneity with in galaxies 
themselves (|Niindl201lt ICampisi et al.l)2009D . Still, the 
apparent amplification in the GRB-to-SFR rate in very 
low-mass galaxies relative to moderate-mass galaxies is 
more suggestive of a gradual trend, and the direct detec- 
tion of a handful of GRBs in high- metallicity enviro n- 
ments (|Pr ochaska et al.l I2009L iLevesque et al.l l2010bl lal: 
iKriihler et al.ll2012al ) is difficult to explain with a sim- 
ple cutoff. 

The apparently improved consistency of the GRB 
and SFR toward z ss 2 could reflect an overall de- 
crease of metallicity in luminous galax ies compared to 
similar-mass objec ts at lower redshifts ()Erb et al.l 120061 : 
iReddy et al.l [2010). The speed of this transition is 
notable — GRBs in the quartile of star formation repre- 
senting the most massive galaxies are common at z = 
1.5-2.0 but almost entirely absent at z < 1.5. This 
rapid trend hints at a sharp change in the properties 
within the Universe's most massive galaxies at this red- 
shift. This could be interpreted as a consequence of the 
gradual buildup of metals in the most massive galaxies 
pushing the metallicity of these systems above a criti- 
cal value beyond which GRB production is greatly sur- 
pressed. Alternatively, perhaps this epoch in cosmic his- 
tory corresponds to the end of the accretion of streams of 
metal-poor gas onto massive galaxies (permitting chem- 
ical homogenization and a sharp drop in the amount of 
metal-poor star- formation) . 

Independent of mass, GRBs may exhibit a secondary 
preference for systems with high sSFRs beyond that 
predicted from their higher overall SFRs. This could 
also be an imprint of metallicity: hig h-sSFR galaxies 
have lower metallicity at a given mass ([Mannucci et al.l 
1201 lh . although the surprising strength of this effect in 
the highest-sSFR system requires (at minimum) a very 
strong metallicity dependence to explain, and may also 
require consideration of additional factors affecting the 
GRB rate other than metallicity. 

As an alternative to metallicity, large variations at 
the top end of the IMF in different galaxies could also 
produce significant variation in the GRB rate by affect- 
ing the relative numbers of extremely massive stars that 
(presumably) produce GRBs relative to the smaller (but 
still massive) stars which dominate the long-term UV 
luminosity. However, if indeed the ratio of GRB rate 
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to SFR varies across the entire mass range (and given 
the need to produce an order-of-magnitude change in 
the rate of production of the most massive stars), this 
effect would have to be so pervasive — or limited only to 
the most massive stars — that it is difficult to imagine 
such an effect would escape notice by other means. A 
more subtle effect to recognize would be a change in the 
binarity properties of the initial stellar population, an 
effect which would be observationally difficult to recog- 
nize directly but which could affect the production of 
GRBs dramatically if the progenitor is a binary system. 
Nevertheless, the simplest interpretation of our results at 
present is a direct dependence of the formation rate of 
the progenitor on metallicity. 

Our conclusion — a GRB rate (relative to that of star 
formation) that is a strong function of host-galaxy 
environment — differs from that of some other recent work 
reexamining the GRB rate in the context of the mass- 
sSFR-metallicity relation (jMannucci et all 1201 lh . We 
note that these authors considered only the explanation 
of the apparent observed metallicity bias itself and did 
not address in detail the observed mass distribution of 
GRB host galaxies beyond comparison of a single mass 
bin to a single study. While it is natural to expect an 
SFR-selected population to skew toward lower masses 
and to higher (s)SFRs (and therefore metallicities) than 
a galaxy-count-selected populati on, the magnitude of the 
effect should be relatively small (jKocevski fc Westll2011f) 
and has difficulty explaining the frequency of GRBs in 
extremely low-mass en vironments (and the low rate in L* 
galaxies.) In addition, iGraham fc Fruchterl ()2012f ) have 
recently considered predictions for the mass and metallic- 
ity distributions of GRBs for SFR-weighted galaxy pop- 
ulations and SN hosts, and they continue to find large 
inconsistencies between the expected and observed GRB 
rates in the same z < 1 population unless metallicity 
dependence (or analogous effects) are considered. 

It is also important to note that our result does 
not directly contradict previous work showing gen- 
eral consistency between the GRB r ate and SFR in 
high- redshift host populations (e.g.. iJakobsson et al 
2005; Fvnbo et al.ll2008t IChen et alj feoOQ; Sa vaglio et al 



20091 ). These studies considered only (or primarily) op- 
tically bright bursts and only UV-based estimates of the 
SFR density. The hosts in which GRBs are most dramat- 
ically under abundant tend to be heavily dust obscured, 
so consideration of the obscured GRB population (and 
dust-unbiased estimates of galaxy SFR) is essential to see 
these effects clearly. In addition, since the effect is both 
intrinsically smaller and observationally more difficult to 
recognize at z > 1.5 relative to at z w 1 (and may even 
disappear completely at higher redshifts), it is no sur- 
prise that work focused on higher- redshift populations 
similarly did not show the same trends. 

The clear dependence of the GRB rate on host-galaxy 
properties out to at least z ~ 1 provides a strong cau- 
tionary note about the use of GRBs as direct tracers 
of the cosmic SFRD before the intrinsic reason for this 
variation is better understood. While it is encouraging 
that the magnitude of the trend for GRBs to prefer low- 
mass galaxies decreases at z ~ 2 (suggesting that GRBs 
become better tracers of star-formation as the average 
cosmic metallicity drops), we cannot yet confirm that it 
becomes insignificant at any redshift. One of the most 



exciting broader applications of GRB studies is to probe 
the contribution to the star-formation rate from low- 
to moderate-lumino sity star-forming galaxies at z > 5 
(jTanvir et al J 12012). a population that no other obser- 
vational techniques can currently reach. However, such 
work needs to be supported by better multiwavelength 
studies of GRB hosts at z = 2-3 to more conclusively 
evaluate whether or not the tools we wish to apply to 
these earliest epochs of cosmic history work in the way 
that we expect at more familiar redshifts. 

Even if GRBs do not represent ideal population probes 
of the star-formation properties of the Universe, they 
certainly remain useful as individual probes of specific 
galaxies and their internal structure. The clear correla- 
tion between afterglow-sightline and host-averaged dust 
attenuation indicates that the dust within high-z galaxies 
is usually distributed fairly homogeneously, suggesting 
that empirical reddening corrections used are generally 
reasonable. At the same time there are significant num- 
bers of exceptions (of both obscured GRBs in blue hosts 
and unreddened GRBs in obscured hosts), indicating 
that dust in dense heterogeneous clouds is also present. 
This heterogeneity is most pronounced in intermediate- 
mass galaxies, whereas low-mass galaxies are almost al- 
ways homogeneously dust free and high-mass galaxies are 
ubiquitously dusty (although some sightlines are surely 
much dustier than others). 

On a similar note, we have now determined unambigu- 
ously that GRBs can in fact form in several environments 
that heretofore not been found to host these events: in 
ULIRGs and in very massive galaxies, even if the num- 
bers of GRBs within such systems is lower than might 
be expected for pure SFR tracers (especially at z < 1.5). 
While this study shows that GRB afterglows in these en- 
vironments are almost always heavily obscured, for the 
most luminous bursts in the Universe this is not always 
a limitation, especially if large-telescope spectrographs 
can observe the event rapidly. This has in fact occurred 
for two events already (GRBs 080607 and 070802), and 
while these events are rare, future rapid follow-up ob- 
servations of GRB afterglows should identify more such 
targets, giving us perhaps the only means of probing deep 
inside high-redshift LIRGs and ULIRGs at this level of 
detail. 

In addition, the tendency of GRBs to prefer faint 
galaxies actually has some advantages, since this is ex- 
actly the population that is difficult to understand from 
field-survey techniques. For example, perhaps one of the 
most important results from our study is the fact that we 
did not find any dust-obscured GRBs within very faint 
or low-mass galaxies. A long-standing concern in current 
estimates of the cosmic SFR density is that the presence 
of additional, unaccounted, heavily obscured star forma- 
tion within very faint galaxies cannot be directly obser- 
vationally ruled out. If low-mass galaxies at z « 2 were 
somehow able to conceal a majority of their star forma- 
tion within optically thick regions, this would have pro- 
found implications for the sites of cosmic star formation. 
However, if this sort of star formation were cosmologi- 
cally relevant, we would observe significant numbers of 
highly obscured GRBs from low-mass galaxies. We do 
not see such a population; GRBs in low-mass, NIR-faint 
galaxies are unobscured with only occasional exceptions. 

In addition, if the reason for the GRB rate dependence 
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can be understood (for example, if it can be firmly pinned 
on metallicity), then the variations we see can themselves 
be used as tracers of cosmic evolution. For example, the 
apparent commonality of GRBs in massive galaxies at 
z w 2 compared to their relative rarity at z w 1 probably 
says something about the evolution of the characteristics 
of systems with redshift: massive z w 1 starbursts seem 
to have a higher metallicity (or other systematically dif- 
ferent feature relevant to the GRB production rate) than 
their counterparts at higher redshifts. 

In short, while the apparent noncorrespondence of the 
GRB host population with the simplest expectations is 
disappointing from the perspective of using them as easy 
probes of evolution of the SFR density, further study 
of these intriguing objects nevertheless may continue to 
provide important clues about the nature of the high- 
rcdshift universe. 
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GRB Afterglow Properties 
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( a ) Redshift of the GRB or host galaxy. Host-galaxy photometric redshifts are given to one significant decimal 
place (see Table [4] for uncertainties). 

< b > Gamma-ray fluence in the Swift BAT 15-150 keV band, with the exception of GRB 051022 for which the 
HETE WXM/FREGATE 30-400 keV band is used. 

( c ) Excess X-ray absorption above the Galactic value, expressed as an equivalent neutral-hydrogen column at 

2 = 0. 

( d ) X-ray spectral slope assumed in the afterglow extinction-fitting analysis. The left column is the measured 
value with uncertainties (from the automated analysis of Butler et al.); the right column shows the "minimum" 
value assumed in deriving limits on extinction. It is assumed that /3x > 10. 

( c ' Limiting value (upper limit) of the observed optical-to- X-ray index (corrected for Galactic extinction). The 
optical or NIR point after t > 1000 s with the strongest constraint is chosen. 

( f ' Limiting extinction (Ay in the host frame) based on analysis of individual points relative to an extrapolation 
of the the X-ray flux assuming the spectral index indicated under the fix column. The optical or NIR point after 
t > 1000 s with the strongest constraint is shown. 

( g ' Measured or limiting extinction (Ay in the host frame) based on a fit of at least two optical/NIR points in 
different filters. 
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Afterglow and Host Positions 
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0.067 


04:35:55.03 


-12:40:45.38 


1.0 





.006 


090417B 


13:58:46.590 


+47:01:05.00 


1 


.0 


XRT 


0.017 


13:58:46.65 


+47:01:04.37 







.006 


090709A 


19:19:42.640 


+60:43:39.30 





.1 


P60 


0.090 


19:19:42.71 


+60:43:39.54 


0.9 





Oil 



( a ' 90% confidence radius of the error circle. 

( b ) Telescope providing the most precise afterglow position. See text for references. 

( c ) Foreground extinction, from Schlcgcl ct al. (1998). 



( d ) Minimum aperture radius used for photometry of the host galaxy. 

( e ) Probability of inclusion of an equivalent or brighter galaxy in the error circle due to chance. 



TABLE 3 
Host-Galaxy Photometry 



GRB 


Filter 


Mag. 


Unc. 


Flux 


Unc. 


Telescope 


Reference 


050915A 


9 


= 25.20 


0.17 


= 0.33 


0.06 


Kcck-I/LRIS 






V 


= 25.06 


0.11 


= 0.38 


0.04 


Keck-I/LRIS 






H 


= 24.56 


0.16 


= 0.49 


0.08 


VLT/FORS2 


(a) 




1 


= 23.86 


0.09 


= 0.73 


0.06 


Keck-I/LRIS 






K 


= 20.69 


0.24 


= 3.57 


0.87 


VLT/ISAAC 


(a) 




3.6 


= 18.88 


0.10 


= 7.86 


0.76 


Spitzer-TRAC 






4.5 


= 18.38 


0.10 


= 7.97 


0.77 


Spitzer-IRAC 





(*) Host photometry for all other galaxies will be included in the accompanying 
online supplement. Magnitudes are given in standard calibration system (Vega 
or SDSS) with la uncertainty. Observations without associated references are 
calculated in this work, although in some cases images are taken from previous 
studies — see text for details. Magnitudes are in the Vega or SDSS system and 
not corrected for foreground extinction, while flux values in /iJy are corrected 
for foreground extinct ions . 

( a > iHiorth et all (|2012ft : Malesani et al. 2013, in prep. 
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TABLE 4 
GRB Host-Galaxy Properties 



GRB a 




2 


b 


SFR C 

Moyr- 1 


1O 9 M 




A v c 
mag 


X 2 /dof 


050915A 


2 


.53 




135.8+^-i 
J J w — 48.3 


36.7+} 6 ,-° 

J J ' —10.4 


1 


t-i+0.14 
'° — 0.20 


4.8/4 


051008 


2 


,9o: 


h0.29 
-0.16 


80 4+? 8 , 2 ; 3 

uw — 52.4 


2 2+ 6 -f 

— 1.1 





84+0-? 9 
•°^— 0.34 


5.8/4 


051022 





.81 




26.2+™ 
-6.6 


17.8 + o-n 
— 3.0 





.73+0-1 3 

— 0.17 


18.2/12 


060202 





.79 




^ w — 1.9 


l.l+nf 
—0.1 


1 


.00+°" 

■w — 0.20 


15.5/3 


060306 


1 


.55 




244.0+i? 8 ; 7 
w — 67.1 


8.0+?,? 

— 0.9 


2 


21 +0.14 
—0.12 


11.5/6 


060319 


1 


.17 




0.0+°* 

— 0.0 


6.8+1- 4 

— 1.2 





n7 +0.17 
—0.07 


2.7/1 


060719 


1 


.53 




— 4.0 


13.0+7,°o 5 

— 12.8 





oy+1.10 

"° ' —0.37 


0.5/0 


060814 


1 


.92 




236.7 + ?o'o 

J ^ ' —18.3 


9.8+°! 

^ ^—0.8 


1 


r>o + 0.06 

'-*— 0.04 


3.2/3 


060923A 


3 


nn +0.75 
uu -1.08 


— 33.0 


149.9+ 42 - 3 


1 


56 +0.29 
■' Jyj — 0.27 


10.8/11 


061222A 


2 


.09 




2.7+°-? 

1 —0.2 


+ 11 

-^'^—0.6 





00+°°° 

—0.00 


13.3/6 


070306 


1 


.50 




7.7+2-9 

— 0.9 


— 8.6 





■oo+S-iS 

— 0.00 


18.0/18 


070521 


1 


.70; 


hO.50 
-0.33 


42.2+"'? 

— 10.7 


32.2+iV 

— 4.7 


2 


.19+°°? 

— 0.13 


22.3/11 


070802 


2 


.45 




79.8+,T Q 5 

' ^ w — 44.9 


i.i+°- 7 


1 


■40+0- 34 

• w — 0.30 


3.4/3 


071021 


2 


.45 




107.8+|^ 


98.1+"' 8 


1 


fi o+0.26 
27 


9.7/6 


080207 


2 


.09 




47.2+f, 7 f 

— 44.6 


120.1+I 44 

— 49.3 


2 


,36+o :45 7 

— 0.37 


48.3/15 


080325 


1 


.78 




15.2+f 4 4 


153.4+? 9 

— 25.0 


1 


.09+0? 8 

— 0.13 


18.4/8 


080607 


3 


.04 




18.3+|° 


36.2+|? 


1 


,08+°-i 3 

— U. 11 


7.8/3 


081109 





.98 




51.1+1-1 


8.8±i;i 


1 


oc+0.08 
— U.Uy 


19.3/15 


081221 


2 


.26 




73.0+H1 


65. 9±\\ 


1 


oq+0.09 

• OiJ Q Qg 


29.5/9 


090404 


3 


nn +0.84 
uu -1.75 




98.9+L 7 o 

— )) 1, Ct 


1 


o 9 +i:i6 

22 


10.8/5 


090407 


1 


.45 




16 - 6+ 15 8 2 


20.1+ 2 1 5 - 8 


1 


.7710^ 


8.3/5 


090417B 





.34 




0-5+gl 

— U.o 


3.4+g-| 

— U.O 





o 7 +0.12 
— 0.08 


21.3/15 


090709A 


1 


o n +0.55 
' 8U -0.66 


8.3+|-| 


10.0+2 5 


1 


-40+EJ'i 


16.0/7 


970228 





.69 




0.5+i^ 

— O.l 


0.3+;;- i 

— O.l 





.63+"-" 


11.3/4 


970508 





.83 




— 0.6 


0.2+°-2 

— 0.0 





.84+0-16 

— 0.19 


8.3/3 


970828 





.96 




35.0+i 2 fi 6 

— 7.6 


0.9+°-? 

— O.l 


2 


-,0 + 0.10 
—0.09 


12.3/4 


971214 


3 


.42 




ss.g+JV/ 


7.1+2-6 

— 2.4 


1 


OC-+0.18 
—0.10 


4.3/3 


980613 


1 


.10 




17 - 9+ 73 


0.4+°? 

— O.l 


1 


.02+0-" 

— 0.19 


20.1/2 


980703 





.97 






5.8+° 4 


1 


10+0.07 
—0.06 


23.5/5 


990123 


1 


.60 




108.2+^-6 

— 50.8 


0.7+° :3 

— O.l 


1 


21+0.17 
—0.19 


4.8/4 


990506 


1 


.31 




0-6+n 2 
— O.l 


26. 5+ 7 :' 

— 21.5 





■oo+Y" 7 

— 0.00 


0.0/0 


990705 





.84 




4.4+°i 

—0.4 


113.0+??- , 

— 18.8 





.00+°°° 

— 0.00 


5.2/0 


990712 





.43 




0.0+°° 

— 0.0 


1.6+°? 

— O.l 





.00+°°° 

— 0.00 


13.7/5 


991208 





.71 




1.0+°? 

—0.2 


0.7+°? 

w ' —0.2 





49+°?f 

—0.17 


4.6/3 


000210 





.85 




o.o+°- 3 

— 0.0 


2.1+°- 3 

— 0.3 





— 0.05 


17.0/6 


000418 


1 


12 




52.4+l 3 „ 7 

— 8.8 


0.8+°? 

—0.1 


1 


oq+0.06 
—0.07 


12.8/5 


000911 


1 


.06 




2.7+ 79 q ° 

— 1.9 


1 2+ 13 
' —0.9 





■ oU n' on 
— 0.80 


1.3/3 


000926 


2 


.04 




8-2 + 3 9 9 9 


4.4+f 8 8 





5 o+0.39 
— 0.29 


0.5/2 


010222 


1 


.48 




0.6+°- 6 

— O.l 


0.7+°- 3 

— 0.4 





.05+0- 33 

— 0.05 


11.6/4 


010921 





.45 




2.7+°i 

— 0.4 


4.1+0-5 

— 0.2 





.48+0-" 

— 0. 10 


21.6/11 


011121 





.36 




1.0+°-; 


13.5+ 4 - 4 

— 4.2 





.00+°°° 

— O.00 


2.2/2 


011211 


2 


.14 




7-0+ n 9 n 4 

— 0.0 


0.1+0-3 

— 0.0 





ig+0.70 
—0.00 


8.1/0 


020405 





.69 




11.6+4-2. 

-2.7 


8.8+11 





.82+0-« 


18.5/3 


020813 


1 


.25 




^'^— 0.2 


9.5tl 3 3 8 





-00 +0 , : n 6 

— 0.00 


6.7/1 


020819B 





.41 




5-8+o's 


84.9+ij 





00+o.os 

—0 00 


43.5/6 


020903 





.25 




0+ 00 


5+°- 2 
u -°-o.o 





cm+o:oo 

■^-0.34 


3.2/0 


021004 


2 


.33 




1 4 8+ 3 - 7 
1 ^-°-2.0 


9 S+ 1 - 
^■ S -0.6 





49 + 0.09 
,yi -0.07 


20.1/6 


021211 


1 


.01 




Q q + 4.6 
8.3_ 7 


2 0+ 11 

z - u -i.o 


1 


7s +0.27 
•'-0.09 


2.4/0 


030328 


1 


.52 




1+51.2 
Z0 - i -15.4 


6+ 35 


1 


06+0- 26 
UD -0.29 


0.8/5 


030329 





.17 




2+ - 1 
u - z -o.i 


l+0° 
u - i -o.o 





co+0.15 

•^-o.is 


11.2/11 


030528 





.78 




6 8+ 4 - 5 
°-°-0.8 


9 1+0.9 
Z - L -1.1 





OQ+0.25 
,uu -0.00 


4.2/4 


031203 





.10 




14 1+0.5 


3+0° 

u - c> -o.o 





04+0.02 

'° -0.02 


239.8/1 


040924 





.86 




9+°- 7 


1 7+ 1 - 6 
1 -'-1.2 





OQ+0.20 
uu -0.00 


1.0/6 


041006 





.71 




3+ 10 


8+ 5 - 1 





00+0.53 
• uu _0.00 


0.1/1 



(a) Burst name. Events below the horizontal rule indicate the pre- Swift 
comparison sample. 

( b ) Redshift. Uncertainties indicate a photometric redshift and correspond 
to the 95% confidence range. 

( c ) Dust-corrected star-formation rate, calculated from the UV/optical/IR 
SED fit. Uncertainties on all parameters except redshift arc la. 

W Total stellar mass. 

( e ) Average dust attenuation. 
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